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^^ ■ A simple and purely algebraic construction of super-energy (s-e) tensors for arbitrary 

fields is presented in any dimensions. These tensors have good mathematical and physical 

properties, and they can be used in any theory having as basic arena an n-dimensional 

manifold with a metric of Lorentzian signature. 

In general, the completely timelike component of these s-e tensors has the mathemat- 
ical features of an energy density: they are positive definite and satisfy the dominant 
property, which provides s-e estimates useful for global results and helpful in other mat- 
ters, such as the causal propagation of the fields. Similarly, there appear super-momentum 
vectors with the mathematical properties of s-e flux vectors. 

The classical Bel and Bel-Robinson tensors for the gravitational fields are included in 
our general definition. The energy-momentum and super-energy tensors of physical fields 
are also obtained, and the procedure will be illustrated by writing down these tensors 
explicitly for the cases of scalar, electromagnetic, and Proca fields. Moreover, higher 
order (super) fc -energy tensors are defined and shown to be meaningful and in agreement 
for the different physical fields. In flat spacetimes, they provide infinitely many conserved 
quantities. 

In non-flat spacetimes, the fundamental question of the interchange of s-e quantities 
between different fields is addressed, and answered affirmatively. Conserved s-e currents 
are found for any minimally coupled scalar field whenever there is a Killing vector. Fur- 
thermore, the exchange of gravitational and electromagnetic super-energy is also shown 
by studying the propagation of discontinuities. This seems to open the door for new types 
of conservation physical laws. 

PACS: 04.20. Cv, 04.50. +h, 04.40.Nr, 02.40.Ky 



1 Introduction 



One of the consequences of the geometrization of gravity is the lack of a well-posed def- 
inition of local energy-momentum tensor of the gravitational field. This is sensible from 
the physical point of view because the geometrization follows naturally from Einstein's 
Principle of Equivalence [74], which is an essential ingredient in most theories that incor- 
porate the gravitational field. But due to the Principle of Equivalence, one can always 
choose a reference system along any timelike curve such that the local gravitational field 
vanishes on the curve, and thus the gravitational energy density must also vanish there. 



This is usually referred to as the non-localizability of the gravitational energy [74]. There- 
fore, any theory of the gravitational field in accordance with the Pinciple of Equivalence 
cannot include the concept of gravitational energy-momentum tensor. 

Nevertheless, there are local tensors describing the strength of the gravitational field. 
For instance, in the prominent theory of General Relativity, Bel M, y, 0] and independently 
Robinson (see for instance the recent Ref.p5|) constructed a four- index tensor for the 
gravitational field in vacuum by exploiting the mathematical analogy of this field with 
the electromagnetic one. The properties of the now famous Bel-Robinson tensor are 
similar to that of traditional energy-momentum tensors: it possesses a positive-definite 
timelike component and a causal 'momentum' vector; its divergence vanishes (in vacuum); 
the tensor is zero if and only if the curvature of the spacetime vanishes, and some others. 
The appearance of the Bel-Robinson tensor seemed more an embarrassment than an 
achievement, though, and has remained as a mistery over the years. It ressembles energy- 
momentum tensors, yet it is not such a thing (it cannot be, as explained in the first 
paragraph). It has four indices, instead of the usual pair. It looks related somehow to the 
energy-momentum properties of the the gravitational field, but its physical dimensions are 
not those of an energy density [Q, |6|, ^, [L9|, 84 ] , and so on []. Actually, as first proved by 
Horowitz and Schmidt |54|] , the Bel- Robinson s-e appears at the first relevant order in the 
definition of gravitational quasi-local masses in vacuum. This has been later confirmed 



many times for the various definitions of quasi-local masses, see 11 , 15, 23, 22, 34, 35, M , 
57|, |§,|93|, Hf. 

In spite of all the interpretational problems of the Bel- Robinson tensor, it is a fact 
that it has been succesfully used in many different applications^], see for a selection [Q, 
JH, HI, H, ||, ||, 1^, H, H |0|, H, 0, ^1- Of at least similar importance is that the 



Bel-Robinson tensor arises as a relevant tool in many mathematical formalisms involving 
the gravitational field. Outstanding cases are: the hyperbolic formulations [39], 33[ of 
the Einstein field equations [|], [I?]], the causal propagation of gravity jl], |l(:7Pl; the 
existence of global solutions of the Cauchy problem [25, 53, 82], and the study of the global 



stability of spacetimes [26, 82]. Furthermore, the Bel- Robinson tensor (or some of its 
generalizations) appears as a pertinent object deserving attention in many other theories 
of the gravitational field, such as supergravity, string theory or the now fashionable M- 



theory [51], as has been demonstrated explicitly several times, see e.g [29, [3C], 31, 32, 



and references therein. All in all, it seems inevitable to ask the question of whether or 
not the Bel-Robinson tensor and its applications can be translated to more general fields 
and theories. 

Generalizations of the Bel- Robinson tensor to non-vacuum cases in General Relativity 



^ue to all these reasons, Bel j|, g, ||, [7| introduced the name super-energy to refer to the 'enigmatic' 
quantities described by his tensor, but wishing to keep some grammatical link with 'energy'. 
2 Deser has applied the epithets ubiquitous [E9| and immortal [B0| to the Bel-Robinson tensor. 



and to other physical fields have been seeked several times. To start with, Bel himself 
provided the generalization to a non- vacuum spacetime in General Relativity [^, || . This 
Bel tensor keeps most of the previous good properties, but it has non-zero divergence in 
general ]|, ^, [l9|, |89], [Xj. This is natural if one does believe in the physical reality of 
the s-e concept, because the existence of matter currents in the spacetime implies that 
other types of super-energies are present. Despite this, some other generalizations trying 
to keep the divergence-free property were built [|(], |28|, [|, but they lacked all other good 
properties of the Bel and Bel-Robinson tensors, and thus they have been hardly useful. 
In fact, as I shall discuss in section [7|, not even their divergence-free property has physical 
relevance, because it does not allow to construct any conserved current or quantity. A 



more fruitful line was launched by Chevreton [ 24 ] , who constructed a s-e four-index tensor 
for the electromagnetic field, showing in particular that in Special Relativity it had all 
the good properties of the Bel-Robinson tensor, including the divergence-free one. The 
importance of Chevreton's idea is that he tried to carry the s-e concepts to other fields (see 
also [5g] ) ; very recently this has also been done for the scalar field in Special Relativity by 



Teyssandier [95] with interesting results. Yet another type of Bel- Robinson generalization 
was put forward by Oktem 75], and a different approach has been pursued by Garecki 
[42, 43, 44, 45, [S], |47j, see [48| for a recent compilation, although their results were only 



partly useful. Recently, Robinson himself [85] has given some clues of how to generalize 
the Bel tensor. Nevertheless, all these results are only of partial relevance and of limited 
application, and have not answered the question of how to define s-e tensors in complete 
generality with the required mathematical and physical properties. 

One of the main purposes of this paper is to prove that the concept of super-energy 
tensor is not exclusive of the gravitational field and that it can be associated in a unique 
and precise way to any physical field. To that end, I put forward what seems to be the 
correct and universal generalization of the traditional s-e tensors. The only pre-requisite is 
an n-dimensional differentiable manifold with a metric of Lorentzian signature, and thus 
it is valid for almost all known theories (including or not the gravitational field), such 
as supergravity, string theory, M-theory, Kaluza-Klein and gauge theories, higher order 
lagrangians, Brans-Dicke and scalar-tensor theories, dilaton gravity, General and Special 
Relativity, and many others. First of all, I give a fully general algorithm to decompose 
any given tensor field into its 'electric' and 'magnetic' parts with respect to any observer 
(section |2|). Then, I present a new and general definition for the basic s-e tensor of any 
tensor field (section ||, Definition |3.1| ) . This definition is very simple and purely algebraic 
and keeps all typical algebraic good properties of the Bel-Robinson tensor. In particular, 
the completely timelike component can be seen as the timelike component of a s-e flux 
vector which is causal, and it is positive definite, as explicitly shown in section y. In fact, 
this component is shown to be the sum of the squares of all the electric and magnetic 
parts of the tensor. 

A specially remarkable property of the Bel- Robinson tensor that has played an impor- 



tant role in some of the above-mentioned applications ]12|, 16, 17, 18 



9TJ is that 



a positive (resp. non-negative) number is obtained when contracting the Bel-Robinson ten- 
sor with any four future-pointing timelike (resp. causal) vectors [80, 18, 12, 39, 9C[. This 
property will be called the dominant property [1^, |8^, |9(| in analogy with the dominant 



energy condition for the energy- momentum tensor [53|. In my opinion, this is one of the 
important requisites to be kept by any good definition of s-e tensor, even more important 
than the formulas related to the divergence of the tensor. One of the main results in this 



paper is the proof that Definition 3.1 achieves this goal, for the dominant s-e property is 



shown to hold in full generality in section ||. This is very important because the cited 
results and applications can be then translated to general fields by using the s-e tensors 
of Definition 3.1. As a relevant example, a recent theorem [pi] in which the dominant 



property of the s-e tensors is essential and that proves the causal propagation of the 
gravitational and many other fields in general spacetimes is included here. 

All the previous properties are also shown to hold for a more general family of s- 
e tensors obtained by permutting indices in the basic s-e tensor. This is explained in 
section ||. However, as I will argue in this paper, only the completely symmetric part of 
these tensors (which is unique) plays a role in the physical applications. This can be seen 
as indirect proof of the uniqueness of the s-e tensors herein presented. 

The physical applications are left for the next sections. In section n, the Definition 



37T| is used to build up and study in some depth the s-e tensors of the gravitational, 
electromagnetic, Proca, and scalar (with or without mass) fields in general theories. I 
shall show that there exists, for any field, a hierarchy of infinitely many '(super) fc -energy' 
tensors involving the field and its derivatives up to a given order related to k. The first 
level in the list depends on the physical characteristics of the particular field so that, for 
typical fields the energy-momentum tensor arises naturally as a first step in the hierarchy, 
but for the gravitational field the first natural step is at the super-energy level. This 
hierarchy can be put in correspondance for the different fields, and moreover it provides 
an infinite set of conserved quantities in flat spacetimes (such as in Special Relativity). 

Finally, the fundamental question of whether or not the s-e quantities can be trans- 
ferred between different physical fields (such as energy does) is considered. The results 
obtained are very promising and illuminating. In particular, I prove that one can always 
construct divergence-free currents in any spacetime with a minimally coupled scalar field 
and at least one Killing vector (a symmetry). These are mixed currents, in the sense that 
they involve the sum of the s-e tensors for the gravitational and the scalar field, so that 
they prove the exchange of s-e quantities between these two fields |)(|. Furthermore, these 
currents reduce to the corresponding generalized Bel-Robinson conserved currents in the 
absence of the scalar field, and to the divergence- free currents of the scalar field previously 
built in flat spacetimes if the gravitational field is switched off. Another enlightening re- 
sult is obtained by studying the propagation of mixed electromagnetic-gravitational wave 
fronts, where appropriate (super) ^-energy mixed quantities are shown to be conserved 
along the wave-front if there are conformal symmetries. These conserved quantities nec- 
essarily involve the gravitational and the electromagnetic contributions, and thus they 
are also a proof of the s-e exchange and its validity at the various hierarchical levels. All 
this opens the door to a new type of conservation physical laws which deserve further 
investigation in the different available physical-geometrical theories. 

2 Tensors as r-fold forms. Electric-magnetic parts 

Let V n be any differentiable n-dimensional manifold endowed with a metric g of Lorentzian 
signature (-,+,. . . ,+). Indices in V n run from to n — 1 and are denoted by Greek small 
letters. Thus, the components of g in any given basis are written as g a p. "Spatial" 
indices will also be needed, so that Latin small letters i,j,... will run from 1 to n — 1. 
Similarly, capital letters I,J,... stand for indices running from 2 to n — 1. Boldface letters 
and arrowed symbols indicate exterior p-forms and vectors, respectively. The covariant 
derivative associated to g is denoted by V and its Riemann and Ricci tensors by R a /3\^ and 



Rp U = R p p pu , respectively (sign conventions as in pq]). The canonical volume element 
n-form is denoted by rj, with components rj ai ... an . Round and square brackets embracing 
any number of indices will denote the usual symmetrization and antisymmetrization, 
respectively, so that for instance n ar ... an = ?7[ a;i ... a „], Rpu = R(^v)- Equalities by definition 
are denoted by =, and the end of a proof is signalled by ■ 

A basic operation in what follows is the standard Hodge dual, denoted by *, and defined 
as follows. Let A™ w ...^ = ^ n hu 1 ...u ] denote any tensor with an arbitrary number of 
indices schematically denoted by {&}, plus a set of p < n completely antisymmetrical 
indices H\ . . . fi p . Then, the dual of A with respect to \i\ . . . /j, p is defined and denoted by 

A m , _ J. A mn...to (1) 

Hp+1-..IM* p\ "^ pn 

where the * is placed over the indices onto which the dual operation acts. Observe that 
A^s * = A^ * . It is easy to prove that 

A fsi} ** =(-ir {n - p)+l A^ fll . 

which can be written simply as ** = (— l s )P( n ~P)+ 1 . Consider now any tensor A^ * ^ ^^ Mi 
which is completely antisymmetric in the sets n± . . . n p and v\ . . . v q and also has any 
number of other indices signalled by {fi}. Then one can form the duals with respect to 
Hi . . . (j,p, or with respect to i/i . . . u q , or to both. The fundamental identity concerning the 
dual operation is given by 



A{n}[vi...v p sVp+l-Vn] _ q ( 2 \ 

si u 1 ...u q ( J^ p+ i...^ n — U, \4) 




* u "+*- Un + 1 ( n ) A^-^u u 8 Vq+1 " Vn , =0, (3) 

where ^\'"^ denotes the Kronecker symbol of order s, see Appendix. The above relation 
together with many others derived from it which will be used later are proven in the 
Appendix. 

Let us now show how to decompose any given tensor into 'electric' and 'magnetic' 
parts. Take any m-covariant tensor L n , and consider the set of indices which are 
antisymmetric with /ii, plus [i\ itself. This set constitutes a block of (say) ri\ < n 
completely antisymmetric indices, which will be denoted by [m]. Take then the next 
index in tu x ...u m which is not in [ni] and do the same, forming thus a second block [7^2] 
of ri2 < n completely antisymmetric indices. Continue like this until all indices of t /il ...^ m 
belong to one of these blocks. Then, we have separated the m indices of t Ml ... Mm into 
(say) r blocks, each containing nx (T = 1, ... ,r) completely antisymmetric indices, with 
n\ + . . . + n r = m. In this way, we can consider t lil ... llm as an r-fold (rai, . . . ,n r )-form 
which will be denoted schematically by £[ m i .r nr ,i where [nr] indicates the T-th block 
with nx antisymmetrical indices. This can always be done because, even if tui...u m has 
no antisymmetries, it can be seen as an m-fold (1,. . . ,l)-form. Several examples of r-fold 
forms are: the volume element ?7 Atl ... A t n is a simple n-form; F^ v = F^y] is a simple 2-form, 
while V pFfiu is a double (l,2)-form; the Riemann tensor is a double symmetrical (2,2)-form 
(the pairs can be interchanged) and the Ricci tensor is a double symmetrical (l,l)-form; 
a tensor such as t^p = tu up -\ is a triple symmetrical (l,l,l)-form, etcetera. 



At this stage, one can define all possible duals by using the * acting on each of the 
[n-f] blocks, obtaining the following tensors (obvious notation): 

t * t t * / * * * 

[n— ni],...,[n r ] [ni],...,[n— n r ] [n— ni],[n— 712],..., [fir] [w— ni],[n— 712], ...,[n— n r ] 

There are 1+1 )+...+ [ =2 r tensors in this set (including £[m],...,[n r ]); each 

of which is an r-fold form (except when ny = n for some T, but these special cases will 
be treated as self-evident). All these tensors allow to construct a canonical "electric- 
magnetic" decomposition of £ Ml ...^ m relative to any observer. To see it, take any timelike 
unit vector u, u^u^ = —1, and contract all the above duals with r copies of u, each 
one with the first index of each block. Whenever u is contracted with a 'starred' block 
[n — ny], we obtain a 'magnetic part' in that block, and an 'electric part' otherwise. Thus, 
the electric-magnetic parts are r-fold forms[j which can be denoted by 

{g EE „ . Jp [ni-l],[n 2 -l] [nr-1]) (it H P -^; E) \n-m-l]<\m-i\ K-lb • • • 

r r-1 

••• ! (u<E ■■■fi^)[nil [rv-i-lUn-ry-ll 1 {ti HH P ■ ■ ■ ^D fn-m-lUn-ns-ll IW-lb ••• > 

r-1 r-2 

(^ E . . . E y HH)i m _i] [n-n r -i-l],[n-n r -l] > ••• > {u flH . . . ff )[ n - ni -i] ; [ w - na -i] [n-n,.-lb 



\u< EE „ . 3 J/J2— Mn] ,"2— fn 2 i— »P2— Pn r — ^lfe-Mn, ,^1^2— ^n 2 .— >P1P2— Pn r U U ■■■ u ; 

r 

^ v x * A i n^ + lA i n^+2---/*Ti)t / lt / 2---t / ri2 i'"iPiP2---Pn r 

r-1 

and so on. Here, £r m i .., [ nr ] denotes the tensor obtained from £r ni i .r nr i by permutting 
the indices such that the first m indices are those precisely in the block [m], the next ri2 
indices are those in the block [712], and so on. There are 2 r E-H parts, they are spatial 
relative to u in the sense that they are orthogonal to u in any index, and all of them 
determine t lil ,.. IMn uniquely and completely. Besides, t t i 1 ... fim vanishes iff all its E-H parts 
do. The simplest way to see this is to take any orthonormal basis {e M } with e*o = u 
and decompose all the components of t\ n A\ nr \ in sets according to whether they have a 
'zero' or not in each of the blocks (notice that blocks cannot have any repeated index). 
This provides 2 r sets which altogether comprise obviously all the components of £ui...u m - 
Observe now that each of these sets coincide with one of the E-H parts: blocks with one 
zero correspond to 'E' and those without zero to 'H'. 

Of course, this 'E-H' decomposition agrees with the classical one of any electromagnetic 
field Ffj, v and its generalizations to the Riemann and Weyl tensors H ||, [?], ^37| 70, 72]. To 




illustrate the decomposition, take any simple p-form S^...^ = S^^ 1. In general, this 



11 



has ] independent components. Given any unit timelike u, the electric part of X! 



3 Of course, if one of the blocks has ny = 1, then when concracting with u this block dissapears, so that the 
electric parts of 1-blocks are scalars, and analogously for magnetic parts corresponding to (n — l)-blocks. Then, 
some of the E-H parts may be f-forms with f < r, and this is to be understood. 



with respect to u is the simple (p — l)-form 

V //12— /%> V //12— /%> V / I/IS— Mp] V / fJ,2—fJ.p 

Obviously, (^-E-j has independent components, which correspond to those 

with a 'zero' among the components S^...^, in any orthonormal basis {e^} with e*o = u, 
that is, to Eoia...ip (remember that i,j,... = 1, . . . , n — 1). Similarly, the magnetic part of 
5] with respect to u is the simple (n — p — l)-form 

Now, (f-ff) has ) = I ' J independent components, which correspond 

to those without 'zeros' in the basis above, that is, to S^...^. The sum of the components 
of the electric part plus the components of the magnetic part gives 

Cn — 1 \ I n — 1 \ / n \ 
P-ij + ( p J-(,j 

which is the total number of components of E. The simple p-form 5] can be expressed, 
for any unit timelike u, in terms of its electric and magnetic parts as 



-1)pO-p) 



/ii.../ip — ^"[/iiU- C/ J/i 2 .../ip] ^ / _ \j 'lpi...p„-piJ.i...fi P u yu 11 ] 



P2...p„~ p 



p p+1 p p+2 ...p n - -' "■""" -^-r^-f^; p 

which can be written in a compact form as 

* * 

£ = -« A (*E) + (-l)^-P) [ w A (|if)] , E= -ti A (f I*) - [« A (§23) 

where A denotes the exterior product. These formulae show once again that the E-H 
parts of 5] fully determine E, and that S vanishes if and only if its E-H parts are zero 
with respect to some (and then to all) u. The invariants asociated with S can also be 



expressed in terms of its E-H parts. For instance, from formula (A. 9) in the Appendix 
one gets 

1 1 1 

— Y W~/"f — (5fn (^ U\P-p + 2-Pn (Sen (•Ep\p 2 ---Pp 

.^PL.-Pp^ - , _ 1 y{u I1 )pp+2-Pn\u 11 ) (fl-l)!* ' W '"'' ! ' 1 " ' 

(4 ) 
The general decomposition for arbitrary r-fold forms does all this on each [nx]-block. 

For any spatial tensor S^ u ,,, p relative to u (that is, orthogonal to u in every index), 
the following abbreviated notation 

[S] 2 = S^.pS^-P > 



pv...p 



will be used. Notice that this is positive and vanishes if and only if the whole tensor S 

is zero. In particular, \~ EE . . . E] 2 , \~H E . . .E] 2 , etcetera are all non-negative and any 

r r— 1 

of them vanishes iff the corresponding E-H part so does. 



3 Definition of the basic super-energy tensor 

In this section I show how to construct a good s-e tensor for any given tensor iui...u m - To 
that end, let us start by defining the "semi-square" (i[ ni ],...,[n r l x t[n 1 ]....,[n r ]) of any r-fold 
(m, . . . , n r )-form by contracting all indices but one of each block in the product of t with 
itself, that is to say 

/, ,\ / TT \ i 7 P2—pni, (T2—CTnr I r\ 

(t X tjAl/il...A r /ir- = I XI JZ _ 1 \\ I ''Aip2.--Pn 1 ,---,A r cr2...0-„ 7 .' : Ml —,/X r • (Oj 

This is a 2r-covariant tensor. With this at hand, let us put forward the following general 



definition [B9, 9C|: 



Definition 3.1 The basic s-e tensor of t is defined as half of the sum of the 2 r semi- 
squares constructed with t[ ni ],...,[n r ] an d °^ its duals. Explicitly: 



T AlM ,.. Ar/ir {0=2 ( t [n.],,W >< *[«i],~.[nr]) A .. A 



+ 



+ [t * x t * +...+ 

[n-m],...,[n r ] [re-m],...,[n r ]/ A lAfl ...A r /x r 



+ ...+ ( t • Xt * + . . . 

[ni],...,[n-n r ] [m],...,[n-7i r ]/ A 1W ...ArMr 



[n-ni],[n-n 2 ],...,[n r ] [n-m],[n-n 2 ],...,K]/ ^ 1W< ^ rMr 



+ + (t , . xi , . \ . (6) 

\ [n-n 1 ],...,[n-n T ] [n-n 1 },...,[n-n r }/ A 1 ^i...A r /ir J 

As we will see, this tensor coincides, sometimes, with the energy-momentum tensor of some 
physical fields, and thus the generic name 'super-energy' tensor may be inappropriate. 
However, in a majority of cases the tensor will be related to higher order 'energies' (see 
section |6|), so that it is preferable to keep the 's-e' adjective. Concerning the adjective 
'basic', this is due to two main reasons: first, for massive physical fields the physical s-e 
tensor will be a combination of several 'basic' s-e tensors (see subsection |6.3| ); second, 
in general a linear combination of basic s-e tensors permutting the indices with positive 
constants will also be a good mathematical s-e tensor. This point will be treated is section 

Notice that any dual of t\ ni \ ,...,[■«,.] gives rise to the same basic s-e tensor (|6|). Therefore, 
one only needs to consider blocks with at most n/2 indices if n is even, or (n + l)/2 if n 
is odd. This will be implicitly assumed sometimes in what follows. Observe that if there 
is any block [nx] with nx = n, then this block 'disappears' giving rise to '0-blocks' after 
dualizing. 

The properties of (g) are the following. 

Property 3.1 The tensor ffl) is symmetric on each (Xy/j l x)-pair, that is, 

J-\lHl...\ r IJ T l*J = ■*-(\im)...(\ r flr) I J 



Proof: Notice that all the terms in the definition (|6|) can be separated in two groups, 
one of them without duals in the block [nx] and the other with duals in that block. Now, 
due to formula ( |A.6D in the Appendix, the sum of these two groups is symmetric in the 
corresponding (AxZ-(x)-pair. As this can be done for any block the result follows. ■ 

Property 3.2 // the tensor £r ni i ... [ n j is symmetric in the interchange of the block [nx] 
with the block [nx'] (nt = ny/), then the s-e tensor ftQ) is symmetric in the interchange 
of the corresponding (Ax/^x)- and (Ax'/^x') -pairs. 

Proof: Assume that t\ nr ^\ n ,],...= t...,[n r ,],...,[ny],... with n-j = ny. Then it is obvious 
that 

...,[n— nx ],..., [n—n T i],... ...,[n— n T i],...,[n— nx],— ...,[nx],---,[n—n T i],... ...,[n— n T i],...,[nx],--- 

and thus, from the explicit expression @ it follows immediately that 
■*■ ...\xfix— A-f'MT'— i J ...A T //^ T /...A-f/iT-- l J 



Property 3.3 If n is even, then the s-e tensor §Q) is traceless in any (Ax ^x) -pair with 
nx = n/2. 



Proof: Assume that n = 2nr for some T. Then, as in the proof of Property p.l| , we 

can separate all the terms in (^) into two groups, one with [nx]-blocks and the other with 

* * 

[n — nx] = [nx]-blocks. Then, application of formula ( |A.8j ) in the Appendix to the trace in 
the (Ax^x^pair leads immediately to 

T.% T ... {*} = 0, ifn = 2n T 

■ 
One of the good properties of the s-e tensor (||) is that its completely timelike com- 
ponent has the mathematical features of an energy density, which can be also seen as 
the timelike component of a vector with the typical properties of an n-momentum vector. 
These good properties are going to be proved partly now, and the rest will be shown 
in section ||. However, the physical dimensions of this timelike component will not be, 
in general, those of an energy density. Thus, following Bel in his pionering work on s-e 
tensors Q, [5], [g, 0], I shall generically use the term 'super-energy' for these concepts. 

Definition 3.2 The super-energy density of the tensor t relative to the timelike vector u 
is denoted by Wt{u) and defined by 

W t (u) = T Xl ^... Xrllr {t}u Xl u^ . . . u K u^ . (7) 

The basic properties of Wt(u) are 

Property 3.4 For all timelike vectors u, Wt(u) is non-negative. Furthermore, 

{3u such that W t (u) = 0} ^ T AlMl ... ArM? .{t} = «=> i m ... Mm = 0. (8) 



Proof: From (|6|) and (|7j), together with the definition of the E-H parts of t we get 

^EE..^E} 2 + ...+ 



^^{(n^K 



+ •••+117 TTT [^£-_S 2 >0 (9) 

\t=i (ji-nr-ly.J ; J 

from where the first part is obvious and one sees that Wt(u) can vanish for some u if 
and only if all the E-H parts of t vanish. But as we saw in the previous section, all the 
electric-magnetic parts of t vanish iff the tensor t is zero, in which case the whole T{t} 
also vanishes. ■ 

Property 3.5 The s-e density is half the sum of the squares of all the components oft 
in any orthonormal basis {e^} with e*o = u: 

1 n— 1 



W t (e ) = r ... {i} = - £ |t w 

pi,...,p m =0 



Proof: It is enough to remember that the E-H decompostion is simply a decomposition 
according to whether a block has a zero or not in any orthonormal basis (see previous 
section), and then use formula (||). ■ 

In fact, there is another interesting formula for the super-energy density which is very 
useful in some calculations, and allows to compare with other similar positive quantities. 
Given any unit timelike u, let us define the 'positive-definite metric' relative to u as 

/V 0") = 9fiu + 2Ufj,u v . (10) 

As this is positive definite, it allows to construct positive quantities when contracting with 
tensors. Actually, we have 

Property 3.6 The s-e density of the tensor t relative to u is proportional to the 'square' 
of t with respect to the positive- definite metric h(u) according to 

W t (n\ — -\T\ If t h^ lUl h^l^l h PlC71 hPnrVnr 

/ \x=i nT 7 
Proof: The case of simple p-forms XI is proved in Property 13.81 below. The general case 



is left for Section M, see Corollary 4.2 



The s-e density can be seen as the timelike component of some 'momentum' vectors 
constructed from the s-e tensor. The precise definition is 

Definition 3.3 The super-energy flux vectors of the tensor t relative to the timelike vector 
u are denoted by Pt{u) and defined by 

rp» (£) = - r Alw ... AT _ lMY _ 1 V..A^{t}w Al « W • • • U A T-^-i^r _ _ _ u A rnMr _ (n) 
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There are, in general, r different s-e flux vectors, one for each (Ax^iO-pah- However, 
sometimes some of these vectors may coincide with each other, such as for instance when 



the Property 3.2 holds for some pairs, in which case the corresponding s-e flux vectors are 



equal. The s-e flux vectors can be decomposed with respect to a unit u into their timelike 
component and the corresponding spatial part as 



: P? (u) = W t (u) u u + (5 + u v u p ) r P t 



T p" <u) 



Property 3.7 Pt (u) are causal vectors with the same time orientation than u. 

Proof: That Pt (u) is causal, that is, either timelike or null, is a particular case of a 
more general property of basic s-e tensors to be proven in the next section, see Corollary 



4.5. On the other hand, we have from ( |il|) and Property 3.4 

r P t p (u)u p = -W t (u)<0 

which proves that Pt (u) is future-pointing if u is future-pointing. ■ 

Let us present some s-e tensors of type (|6|) explicitly. Starting with a scalar / we 
simply have 

T{f} = \f. 
For any given simple p-form E^.,.^ = Sr^ 1 ^ i, definition (||) produces 



I- '"J 2\[p-l)\ F HP (n-p-ly. \p p+2 ...pn 



and then using formula (|A.S|) in the Appendix we equivalently obtain 



T v{ S b]} 



1 



(p-1)! 



| I Zj Ap 2 ...Pp Zj A' 



y p2—p P 



1 

2p 



9\v£ 



PlP2—Pp 



y j pip2---p P 



so that 



TV {S M } = T, A {E M } , T% {S M } = (p - - ) -E piP2 ,„ Pp E^-^ 



(12) 



(13) 



(14) 



These explicit expressions allow to prove Property 3.6 for the case of simple p- forms. 
Property 3.8 The s-e density of the simple p-form S relative to u is proportional to the 



'square' ofS with respect to the positive- definite metric (It), that is 

1 



W E (u) 



9 I ' L -'Pi---Pp- L -'is 1 ...i/ p 



h^ 1Vl . . . hf* pVp . 



Proof: From (|l3|) and the definition of (5-E) it follows 



W^{u) 



(p-1) 



x P2---pp 1 

2p~ 



(H^r*-^-* s " 



P2...pp 
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while using ( |lO[ ) together with the properties that (f-EJ is spatial relative to u and 
completely antisymmetric, one gets 



J PlP2---Pp 



\ / Uo... Un V / 



P2---Pp 



/jM2^2 _ _ _ /jPp^p 



/i 2 •••Mp 
J PlP2---Pp' 



1/3... i/ p 



1/2. ..V p _ 
/jM3^3 _ _ _ /jMp^p 



\ / p 2 ...p p \ / 



which comparing with the previous equation proves the formula. 

Observe that introducing (0) into the formula just obtained we can write 



y y u^ lUl h^v v v 

, £ -'u\...u m *-'v\...v v i 1 ' ...ib 



1 



V 



I ^ in... Up t-> v\... Up 



(n-p-l)\ Ku J ^ +2 -' M 



(%H) (5 m^p+2-Mn 1 

UT"" )Pp+2---Pn\U- rL 1 ' 



+ 



(u E )p2-Pp(u E ) 



fa- l\\\U^JP2-fip 



E rp\p2---Pp 



(15) 



which will be useful in the proof of the general case, see Corollary 4.2 



Consider now any double (p,q)-form if[p],[g] an d take its corresponding if[p],[g] with 



ordered indices: if, 



p 1 ...p p ,v 1 ...v q ^[p 1 ...pp\\ui...v q \- 



Its s-e tensor (p) reads 



T ^ l/fr 1 r il — - I ft" \ f(p'p2-Pp,p02~<Tqt 

-La \ \-"-[p],[g]/ — 2 I ( ]\|/ _ ]_\| - n -ap2...Pp,A<T 2 ...<Tq-"' "T 

J£0P2—Pp,LUTq+2—<7n _i_ 

t 

_| jf J£0Pp+2—Pn,Ll<r2—<7q _|_ 

+ 



1 

(p - l)!(n - g — 1)! ap 2 -p P Mq+2-<rr, 
1 



a 
if 



(U — P — l)l(q — 1)! app+2...pn,A<T 2 ...0- g 

1 



^ t J^Ppp+2—pn,LUTq+2—Crn 

(n — p — l)!(n — q» — 1)! «p p+ *2...pn,Ao- 9+ 2...o-„ 
which after use of formulae ( A.10| - [A~T2| ) in the Appendix is equivalent to 



(16) 



_|_Kf jr p2-pp, <?2-crq j> 

1 - fv ap2...pp,po-2...o- ( j- fv /3 A yap**-, 

n . K j> a P2-Pp,<n<T2-Oq J „ _ £^ j>pip2-Pp,<Tl<J2-<Tq 

i/A/i- n -ap2---Pp,o-i(T2...o"g- n -/3 T </Q/3i/Ap- n -piP2-..Pp,o-io-2...o- g -'»- 



Lctp 2 ...pp,(7io-2...a-g 

and thus in general 



,. , JS-pip2—pp, cr 2 ...a q 

Xf3 ly p 1 p2...p p , \<T2...(T q ly p 



(17) 



Ta^A^ {%],[<,]} = T (a/3)(A M ) {•*%[?]} • 



(18) 



One can find similar explicit formulas for the basic s-e tensor T a pxpTv \ -Afpiufs] \ of 
general triple (p, q, s)-forms in a straightforward way, and so on. For later use, the form 
of the s-e tensor for a triple (p, q, s)-form Ayrju is presented: 



TapXp.ru { A \p],[q],[s]} 



1 



(p-l)!( 9 -l)!( fl -l)! 
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A, 



ap 2 ■ . .p p , A<T 2 . . .ff q ,tC,2 ■ ■ -C,s A 



A 02— pp, <T2—(Tq, C2— Cs 



^^PP2-"Pp,Aa- 2 ...<Tq,T<;2"-Cs Y1 a: P v ' 

_l_A /- /- Aa p2 --- p P'\ (72 --- a 1' C2-G 1 

_|_4 , . . 2„P2...pp, (T2— <T„ C2-C 

„ a f I . , , APl-Pp: o-2-O-q, C2-G1 

</a/3 ^pi...pp,Acr2...cr 9 ,rC2---Cs^ 1 P " ^ 

4 . , , APi—Pp> ff 2— c g , Ca— C» 

^ 1 pi...pp,Ao-2...cr 9 ,l'C2---Cs^ 1 P r 

/,, fl - . / 4, Q P2---Pp,0'l...(Tq, C2---Cs_|_ 

i/Ap ^^ 1 ap2...p p ,o-i...cr 9 ,rC2---Cs^ 1 /3 f ' 

A a a 4 „P2— Pp,o-i...o- 9 , C2— Cs 

^ 1 ap2...Pp,<Tl...CT«j,^C2---Cs^ 1 P T 

n (l n a a Aa P2 - p J" "2---o- 9 ,Cl---Cs_|_ 

Utv \^-ap2...p p ,\(T2--aq,C,i--C,a- t± P P ^ 

A a a Aa p2 - pp ' , CT 2-o- 9 ,Ci---Cs N \ 1 

^ap2...Pp,po-2---o- 9 ,Cl---Cs yl P A jT 

-I o „„, i > - 2m-Pp,<n-<rq, Ca-Cs 1 

»o 1 "' Pp ' cri '" CT ' J ' T 

1 

J n art 4 x , t A pl - p P' ^...(TgXl-Cs 1 

j „, 4 , . 4 „P2-Pp,o-i-o-„Ci--Cs _ 

„ r,n\ a A r /■ 4Pi---Pp' cr i--- cr 9.Ci---C s 

2nas p 9 

It is noteworthy that expressions (|I3|), (fl7|) and (|l9|) do not depend on the dimension 
n explicitly. This is in fact a general property, so that a formula for the s-e tensor (g) 
can be given, in general, without any explicit dependence on n. It is enough to expand 
all the duals. Actually, this general expression can be easily guessed for general r-fold 



(19) 



(n\, . . . ,n r )-forms from the cases given in (13), fll7|), (|19|). However, this general formula 



is cumbersome and adds very little to what has been said, so that it will be omitted here. 

4 The dominant super-energy property 



Properties 3.4 and 3.7 are simple particular cases of a quite general and much more 
important property of the basic s-e tensors @, which I call the dominant super- energy 
property (DSEP). The general definition is 

Definition 4.1 A ranks tensor T fJ-1 ... fJ-g is said to satisfy the dominant super-energy prop- 
erty if 

T^M 1 ...#'>0 (20) 

for any future-pointing causal vectors k\ , . . . , k s . 

The simplest examples of tensors with the DSEP are the past-directed causal vectors and 
their tensor products. Another simple example is minus the metric tensor, and any tensor 
products of it 

-9, 9® 9, -9®9®9, , (-l) k g®...®g . 

" v ' 

k 
Actually, we have 
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Property 4.1 The tensor product of any number of tensors with the DSEP satisfies this 
property too. Furthermore, any linear combination with non-negative coefficients of ten- 
sors satisfying the DSEP has also the DSEP. 

Proof: Trivial. ■ 

Let us remark that, in some particular cases, even more general linear combinations 
(for instance with some negative coefficient) of tensors which satisfy the DSEP may also 
keep this property. Property |4.1| shows that the set of all tensors with the DSEP at any 
point of V n has the algebraic structure of a vector 'semi'-space, given that one can only 
multiply by the scalars of the 'semi '-field IR + . 

The DSEP is a generalization of the so-called dominant energy condition for energy- 
momentum tensors |53). The justification for its name is given in the following 



Lemma 4.1 If a tensor Jj^.. „ satisfies the DSEP, then 

T ...o > |r w ... Ms | V>i,...,/i s = 0, ... ,n- 1 

in any orthonormal basis {e u }. 

Thus, the completely timelike component o£Tu u ..u a (which may be called its 'super-energy' 
relative to e*o) 'dominates' over all other components of T^...^. 

Proof: Let {e u } be any orthonormal basis with a future-pointing e*o and define a 
set of (n — 1) + (n — 1) = In — 2 future-pointing null vectors l i by means of 

f = eo±ei, g(£fj+)=g(ir,i-)=0. 

As initial step, from (^(J) one immediately obtains in the given basis 

T All ...„,e£ 1 ...e(f'>0 <=> T ...o>0. 

Next step is to contract with at most one of the null vectors t i . For instance, from (|2 
one derives 

r M1 ... Ms e^ . . . e^-Hf^ > «=► T0...0 ± r ... w > ^ r ... > |t ... w | , 

T^.,.^4 1 ■ ..if^etf > «=► T ...o ± r „.io > <=► T ... > |T ... i0 |, 
and so on for the components with (s — 1) zeros. Next, one computes 

T c^ 1 c ^ s -' 2 (f+Hs-lptPe 1 p-Pe-l //TVs\ 

- t Mi---^s e o • • • e o \ l i l j 'h l j ) 
which is non-negative due to ([H]) and equivalent in the given basis to 

Tq...q ± To.. ,oj + To...jo ± To...jj + T0...0 T To...oj — To...jo ± lo...^- > -4=>- 
•^=^ T0...0 ± To...ij > -<=>■ T0...0 > |To...jj| 



and similarly for any component with (s — 2) zeros. Next one can calculate 
which leads by (p0|) to 



T0...0 > |T ...fcjj| 
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and analogously for the components with (s — 3) zeros. One proceeds in this manner in 
order, and at the (I + l)-th step one uses the schematic combination given by 

if (combination of Z-th step)^ i + £~ (same combination interchanging ± <-> =p)j i 

plus the necessary number of e*o's to complete the s vectors. This finishes the proof. ■ 

A possible generalization of definition (||) that one is immediately tempted to adopt is 

to allow for different positive weights in front of each of the semi-squares on the righthand 



side. This will obviously keep the Property |3.4|. However, Property 3.7 and the DSEP 



would be lost as follows from the following result, which is also the essential step in the 



proof by induction of the fundamental Theorem 4.1. 
Lemma 4.2 Among all the tensors 



Vfa;*,*} s T^TTT^v,.,,^ "* + T^TW* - £«-"" (21) 



2---Pn 



with arbitrary constants ci,C2, the only ones which satisfy generically the DSEP are pro- 
portional to the basic s-e tensor jjiTW), that is, precisely those with c\ = c 2 > 0. 



Proof: By using eq.( [A.9 ) in the Appendix, a calculation similar to that leading to (13) 
proves that 

7V{£ H ; Cl ,c 2 } = ^±^ls Ap2 ... Pp S/-^ - ^ g x „Z piP2 ... Pp ^-e* . (22) 

Let u be any unit future-pointing timelike vector and let v be any causal vector with the 
same time orientation: u^v^ < 0. Obviously, one can write 

v = au + b, a = -u% > 0, ?/% = 0, b 2 = b^b^ > (23) 

and the condition that v be causal reads 



v^v 



V 



-a 2 + b 2 < ^^ b 2 < a 2 . (24) 



Starting with 



one immediately gets for generic S 

ci > 0, c 2 > . 

Concerning the general result, by using ( |22| ) and (|j) we obtain 



+ -((^wl(^)l 2 -^KHl 



(n-p-l)! L V" / J (p 

(^K^)l 2 + I^fr^i(i«)] 2 + |^v (H" " ^-* 
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In order to simplify the calculation one can choose an orthonormal basis {e M } with u = 3q 
and b = be\ with b > 0. Keeping in mind the notation that Latin small letters i,j, . . . run 
from 1 to n — 1 and capital letters I,J,... run from 2 to n — 1 we get 



r v {E b] ;c 1)C2 }«V = a Cl Yl [$ E ) i ... i ) + 



%2<...<lp 



J2<---<jn-p 



z ^ \\ u Jjo...j n _ r , T>-1)\ a 



]2—]n-p 



(ci + c 2 ) b 



(p-1)! a 



'\h--h \u 



iE 



h-i P 



( c ..e.((H 1 „..J 2 + ^i: / ((^), 2 J 2 + 

\ J 3 <-<«p V P/ I 2 <-<Ip V P/ 



+C2 E 

J3<---<jn-p 



rH 



lJ3-jr, 



) + C2 E 



J2<-<Jn — p 



* H \,.,„J + 



l2,---)Ip z / z J2>---,Jn-p 



+(* + < £ (-D" +l " (I*),,..,, (= ff ) 



l2<--.<Ip, J2<--<Jn-p 

where (— 1)'°" denotes the sign of the permutation 

(h • • • -Zp ^2 • • • J-n-p 

2 . . . . (n-1) 

By re-orginizing in perfect squares we arrive at 



l2---Ip ^ ' J2---Jn— p 



T v {£ [p] ; Cl ,c 2 }uV = a Cl £ f(l S ) H ...J + c 2 E ( (*#) 



«3<...<« p 



J3<— <J?ra-p 



h,---,Ip¥ : J2,---,Jn--p , U / \ 

72<---<Jp,J2<---<Jn-p 



+ 



/ 2 -../p 

(.2 



C2 - 4^ C1 I 1 + ^ J 



E ((!*) 

- J2<-<Jn-p 



lis. "it. 



J 2 ■ ■ "Jn — p 



Jl-.-J-n- 



+ 



Due to (p4|), the factor in square brackets satisfies 



62 f 1 + C2 



> 



Cl / C 2 

c 2 - — 1 + — 

4 V ci 



.fiA_£2 

4 V ci 



and thus 



C2 E 

i3<— <in-p 



^{S W ; C1) C 2 }„V>« Cl £ ((§*) ) 



l2,---,Ip^J2,---,Jn—p 

+ci E 

I-2<...<Ip,Jl<...<Jn-p 



lJ3-jr, 



J 2 ■ ■ ■'Jn — p 



Cl / C2 

4 V ci 



£ ((?*) 



J2<-<Jn-p 



J2-..Jn- 
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From here, keeping in mind that the result must hold for generic S, if c\ ^ C2 the 
last term is strictly negative, and given that all the other terms can vanish for some S 
and that the equality can be achieved for b = a, it follows that T\^ \ Eu; c\,C2 > u x v^ 
is certainly negative for some X! if c\ ^ ci. On the other hand, one immediately sees 
that if c\ = C2 the last term vanishes and the righthand side is a sum of squares, so that 
T\^ \ Sfpi; ci, c\ \ u x v^ > 0. This proves the result except for the minor point that u was 
assumed to be timelike from the beginning. However, as this is valid for arbitrary u, the 
result also follows by continuity for null vectors. ■ 

An interesting corollary of this result is 

Corollary 4.1 For all timelike vectors u and w with the same time orientation and any 
simple p-form 5] one has T\^ I Er^i > u x w^ > 0, that is 

^ (5 £ ) (SB)** + "_ j (|ff) {lH)"« -*• > . 

(J)— 1J! V ' P2-Pp \ ' [n—p—\)\\ /pp+2-.-Pn V / 



Proof: This follows at once from the proof of Lemma |4.2| by using c\ = C2 = 1 together 
with b 2 < a 2 (strict inequality). Then, as long as 5] is not zero, T\^ I £r„i > u x w^ is proved 
to be a sum of squares which cannot vanish simultaneously. ■ 

Notice that, in general, neither (%e) (%e) 2 ' P nor (%h) (%h) 

are positive by themselves, nor any other combination of them with weights different from 
those of the Corollary |4.1|. 



Theorem 4.1 The s-e tensors ([^) satisfy the DSEP. 

For General Relativity (n = 4) , the proof of this result has been given in an elegant manner 



by Bergqvist [14] using spinors. Apparently, the use of Clifford algebra techniques may 
also provide a proof of this theorem [J77[] . 

Proof: The proof is by mathematical induction on the natural number r of the 



r-fold forms. From Lemma 4.2 the result holds for r = 1. Thus, let us make the induction 
hypothesis that the result holds for the basic s-e tensor of (r — l)-fold forms and try to 
prove the result for r-fold forms. Same notation as in Lemma [4.2| will be used so that 



u is any unit future-pointing timelike vector, v is any future-pointing causal vector and 



relations (|23|) and (24) hold. Furthermore, the orthonormal basis {e M } with u = e*o and 
b = be\ with b > is also used as well as the notation for Latin small and capital letters. 
Let ty nt v.„An r \ b e any r-fold (m, • • • ,n r )-form and T\ lfM1 ...\ rfMr {t} its corresponding 
basic s-e tensor @. The idea is to prove that T\ 1/11 _\ r _ lllr _ 1 \ rflr {t}u Xr v flr is a linear 
combination with non- negative coefficients of basic s-e tensors of type @ for some (r — 1)- 
fold forms. To that end, define the following set of (r — l)-fold (m, . . . , n r _i)-forms 

y^ 1 ...fi, n ,...,u 1 ...u„ r _ 1 [i2, ■ ■ ■ ,in r ) = tfj, 1 ... flnit ...^ 1 ..M„ r _ 1 ,p 1 p2...p nr u e i 2 ■■■ e iJ r (25) 

Z u , u u, v (*2, • • • , in-n ) = i . u Pl ef . . . ef n "" r (26) 



p,\...p,n 1 ,--;Vl---Vn r _ 1 ,pip2---P 



n — n r 



which contains I + I = different (r — l)-fold forms. From (| 

\ n r LI \ n r J \ n r J 

we have 



r AlM1 ... ArMr {£}u A 'V" = - (< [ni]j .„ l[np ] x <[«i],.-,[n P ]) A x _ -Ar 



+ 
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+ [t , x t * I + ..- + 

[n-m],...,[nr] [n-m],...,[n r ]J A 1/xl ...A rAtr . 



+ ...+ [t . Xt * I +... 

[ni\,...,[n-n r ] [ni],...,[n-n r }/ ^/^...Ar/Ar 



+ . . . + ( t * , Xt, * I -}- . . . + 

[n-m],[n-n 2 ],...,[n r ] [n-ni],[n-n 2 ],...,[n r ]/ \ lfll _\ rflr 



u Xr v^ r 



+ + 1 1 . * xt * 

[n-ni],...,[n-n r ] [n-m],...,[n-n r \/ \ lfll ...\ rflr 



Let us concentrate on the first term of the righthand side. By using the notation intro- 
duced in (25-26) and definition (||) we get due to (p3|) 



(*[ni]..«,[«r] 



xt 



n 1 \,...,[n r 



n 

Al/Ul—Ar-i/Ltr-iArMr \T— 1 ^ nT ~~ ^ 



u Xr v lXr 



'■'■ 2-*/ ^Al/"2— Pn 1 ,».,A r -l<Ta...O"n r _ 1 V ? 2> • • • j W j i/ii 1 '...,/i r ._i r (,*2> • • • > W J + 

i2,—,in r 



' 7 j * Alp2---Pn 1 viAr-10"2---0", 

«2vi*n r 



r _j V-'ii ■ ■ ■ i l n r ) tfi x ni '...,fj, r _ 1 ''" 1 'fi r i 2 ...i 7lr b 



and taking into account that, in the second summatory, only the values of *2> •••)** 
different from 1 are relevant and that 

(-l)-+H6z w »»^... i/lr .^^«r-i (J 2 ,...,J n _ n J for J 2 ,...,J n _ nr ^I 2 ...In T 



where (—1)'°"' is the sign of the permutation 

( h ... In r h 
[2 . . . 

we arrive at 



• • <Jn—n r 

■ (n-1) 



f[m],...,K] x *[ni],...,Kj , , 



u Xr v^ r 



a E (^ni],...Jn P -i](*2>--->»nr) xY in 1 ],...,[n r _i](*2 ) ...,*n r )), 

. v ' Ai/ii...A r _i/i r _ 1 

12<-. .<«n r 

+ E (-i) n+ H6(y M ,... )K _ l] (/ 2) ...,/ nr )x 



+ 



h<---<Inr< J2<---<Jn-n r 



Z [ni],...,[nr-l] ( J 2,---,Jn-n r )j 



or equivalently 



(*[ni],..., 



n r ] x *[rn],...,[n P ] 



Ai/ii...A r _i^ r _iA r ^ r 
E ( 5/ ini],...,[n r _i] (1) ^3> • • • , 4.J x Y[ni],...,[nr-l] C 1 ' h, ■ ■ ■ , In 

h<...<In r 



u Xr v^ r 



Ai/Ui...A r _i/U r - 



+ 
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+ E 

l2<---<In r , J2<---<Jn-n r 



( ^[ni],...,[n. r _i] (h,--- Jn r ) X 

y[ni],...,[nr-i} (h,---,In r ) + -(-1)™ +|ct| -^[m],...,[n r _i] («^2, • • • , J n -n r ) J 

a J / Ai/xi...A r _i/x r _i 

By doing the same with all the terms and using (|6|) one can write 
T\ lfll ...\ rfM , {t}u Xr v^ r =al Y, T AlMl ... A? ._ lMr _ 1 {Yjm],...,^.!] (1, J 3 , • • -,-fnj} + 

[/ 3 <-<Inr 

/ , ^Ai/Ji—Ar-i/Ur-i |-^[m],...,[n r _i] (1) ^3) ■ • • ; Jn-n r ) J + 

f y[ni],...,[n r _i] (^2,---, V) x 



J3 "C- ■ ■ < ~~~'Jn~-n' l 



E 



l2<-<In r , J2<---<Jn-n r 



Y [nxl...,[nr-i\ {h, ■■■Jn r ) + ~( _1 )™ + Z [ni],...,[n r -i] ( J 2, ■■■ , Jn-n r ) 

+ (y% , (h,...,In r )x 

\ [nij,...,[7»r-lj 

Y- (/ 2 ,...,/ n? .) + -(-ir + Hz (J 2 ,...,J n _^) 

[nij,...,[n r _ij a [ni],...,[n r _i] 

+ + all duals + 

+ ( Z [n 1 },...,[n r - 1 ] (J2, • • • , Jn-n r ) X 
Z[ni],...,K_i] (J 2 , • • • , Jn-rv) + -(-1)™ + ^[m],...,^!] 0*2, •• ■ , InJ 



Z» (J 2 ,... , J„_„ r ) X 

[nij,...,[n r _ij 

*• , i (J 2 ,...,«/n-nJ + -(-ir +H y. r (/2,-,/J 

[nij,...,[n r _ij a [rnj,...,[nr_ij 



+ 



Ai/^i...A r _i/i r 



+ 



Ai/ii...Aj.-1/i,. 



+ 



Ai/ii...A,-_i/i r _i 



+ 



Ai/ii...A r _i/j r _i 



+ . 



+ all duals 



which after reorganizing a little bit becomes finally 



r AlM1 ...A rMr {t}« Ar ^ = a < ^ 2A 1M1 ...A P _ 1A t r _i {yin 1 ],...,[n P _i] (1.^3, • • • , V)} + 

[j 3 <...</„ r 

+ /_^ -L\i(i.i...\ r -i(i r -i |-^[ni],...,[n r _i] (.I) ^3; • • • ! Jn-n r ) J + 

l2,---,In r ^J2,---,Jn-n r 



+ E 

l2<-<In r ,J2<-<Jn 



T 



Aijui...Ar- l/Ur-i 



y(/ 2 ,...,/ nr ) + 



+(-lf+H-Z(J 2 ,...,J ri _ nr 



+ 



Ini ,..., n r _i] 
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+ ( 1 ~~ a? 1 ^ T *i»1:.>*-1I*t-1 {■ Z '[ni],...,[n 7 ._i] (J2, • • • , Jn-n r )j > • 

This is a sum of basic s-e tensors for (r — l)-fold forms with non-negative coefficients (due 
to (p3|-P4|)), and thus by the induction hypothesis it follows 

?Ai^i...A r _i/v_iA r ^ 7 . W kx k 2 ■ ■ ■ ^2r-3^2r-2 U r V r >Q 

for all future pointing causal vectors k\, kt, ■ ■ ■ , &2r-2- As u and v are also arbitrary, the 
theorem follows. ■ 

With this result at hand, the proof of Property [Oj is very simple. 



Corollary 4.2 Property \3.q holds 



Proof: From Property [Oj, the result holds for r = 1. Let us use again the mathematical 
induction on r, and assume that the result is true for (r — l)-fold forms. From the proof 



of Theorem LI , but with v = u (which is equivalent to setting a = 1 and 6 = 0), it follows 
that 

7A im ...A rMr Ww Ar U Mr = Y, TXi^i-Xr-ifM-l { Y [ni],...,[n r -i] (*2. • • • ,»!*)} + 

42 <...<in r 



7 j -L\iin...\r-ifir-i j^[ni],...,[rv„i] \j2i ■ ■ ■ ,Jn-n r )j 



32<---<3n-n r 



so that 



W t (u)= Yl W Y{t2 _ lnr) {u)+ Y W z(j2 _ jn _ nr) (u) 

l2<---<in r J2<---<jn-n r 

and using here the induction hypothesis one can write 

/ , ^pi...p„ 1 ,...,pi...p„ r _ 1 (*2) • • • ) 1>n r ) *v 1 ..M ni ,...,cri...o- nr _ 1 \ l 2: ■ ■ ■ i l n r ) + 

1 «2<---<*n r 



"1" 2-~i ^pi. ..fi ni ,..., pi. . ■pn r _ 1 \J2j ■■■ ijn-rir) ■6vi..m„ 1 ,.. 



■ ,u l_...u n 



_j \32i ■ ■ ■ i3n—n r ) I • 



J2<...<>, 



A calculation similar to that leading to (|l5|) and the use of the definitions (25-26) proves 
then the desired result. ■ 



A combination of the proofs and results of Lemma |4.2| and Theorem |4.1| leads easily 
to the following 

Corollary 4.3 Among all the tensors that one can form similar to fiQ) by adding the 
different terms of the righthand side in ffl) multiplied by arbitrary constants, the only 
ones which satisfy the DSEP are precisely those proportional to the basic s-e tensor ^)M 
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Corollary 4.4 For all future-pointing timelike vectors wi,W2, ■ ■ ■ ,W2 r , the strict inequal- 
ity 

r AlM1 ... Ar _ 1/v _ lArMr {t}« . ..< r _!< > o 

holds as long as £r ni ] ... r nr ,i does not vanish. 
Proof: This is similar to Corollary [4.1| . 



Now, it is very easy to prove the first result of Property 3.7. 



Corollary 4.5 The s-e flux vectors (11) are causal with the same time orientation than 
u. 



Proof: By Corollary 4.4 and definition (|Tl|), Pt (u) satisfy 

r Pf (u) w p <0 

for all possible future-pointing timelike vectors u and w. ■ 

This result seems to indicate a causal propagation of the corresponding super-energy 
if one interprets Pt (u) as flux vectors. However, as have been proved, a s-e tensor with 
the DSEP (and therefore with causal flux s-e vectors) can always be constructed for any 
given field tai...u m - This cannot mean that causal propagation is 'universal', because one 
can very easily construct tensor fields iu 1 ...u m by hand which do not propagate causally. 
What happens here is that the super-energy does propagate causally in the sense that 
T Pt (u) are causal, but the real field i jUl ... Mm does not necessarily propagate causally. In 
order to link these two propagations one needs a further condition which is related to the 
divergence of the s-e tensor (||) and thereby to the derivatives ('field equations') of the 
field t fMl _ MMf j m . For instance, in |l6|] we have been able to prove the following interesting 
theorem. Let C be any closed achronal set in V n and -D(C) its total Cauchy development 



(see []53| |8S|, |97J for definitions and notation). 

Theorem 4.2 If the s-e tensor ^) satisfies the following divergence condition 

V P T^- A ^'' {t}w^ . ..w Xr w„ r < f T x ^- x ^ {t}w XlWfll . ..w Xr w Pr (27) 

where f is a continuous function and w = —dr is any timelike 1-form foliating D(Q) with 
hypersurfaces r = const., then 

t/J,l...flm\( = U ^ tfll-.-flm \d(Q = 



Proof: See p| . 



In fact, condition ( |27| ) can be relaxed substantially by allowing 'powers of T" greater 
than one on the righthand side. This theorem allows to prove the causal propagation of 
massive and massless fields in General Relativity, as well as in higher dimensional theories, 



including that of the gravitational field, see [16> 18 1. Furthermore, the uniqueness of the 



solution to the field equations in -D(C) also follows from theorem 4.2, see |16||. Let us 



remark that a key point in the proof of Theorem |L2] is the DSEP for the s-e tensor, so 
that the universality of the construction of the s-e tensor (^|) for arbitrary fields may lead 
to very general results. 
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5 The general super-energy tensor 

As remarked after the definition of the basic s-e tensor (||), the word 'basic' was used 
because more general tensors constructed by permutting indices in j™iw— * r/ir are also 
good s-e tensors. This freedom will be studied in this section. 
The first thing to note is that every tensor of type 

■*|Ul/*2— M2r- 1^2r i*J = ^ M<r(l)/V(2) ■••A t <r(2r-l)M<r(2r) W 

where cr(l), . . . ,<r(2r) denotes any possible permutation of l,...,2r, satisfies the good 
properties shown so far for @ including the DSEP. Thus, let us make the following 

Definition 5.1 The general s-e tensor oft is defined by 

^l/i 2 ...^2r-1^2r 1>J = /_^ C <J ^ M<r(l) /V(2) • • ■M<r(2r- l)PW(2r) V^°) 

w«£/l non-negative coefficients c a . 

Notice that, in general, due to the Property [O] not all the permutations above will 
give rise to a new tensor. Thus, it is very simple to check that there are only (2r)!/2 r 
meaningful independent constants c a or, in other words, T is a (2r)!/2 r -parameter family 



of s-e tensors for t. If the Property ^2 also holds for some blocks [ny] and [n~f], then 



the number (2r)!/2 r is further reduced. Thus, if there are f < r blocks [rax] which can be 

(2r)! 

interchanged, T reduces to a --parameter family. 

2 r r! 
Concerning the properties of T, let us first note that, in general, it does not have any 

symmetry property, nor any particular traceless property. The only thing that can be 

said is that, if the basic s-e tensor is symmetric with respect to all (Ar^iO-pahs, then the 

(2r)! 

general s-e tensor T keeps this property. In this particular case, T is a = (2r — 1)!!- 

parameter family of s-e tensors. With respect to the general super-energy density Wt('u) 
defined as the total timelike component of T with respect to u, we obviously have 



w t (tZ)= (y2cAw t (u) 



so that the original function Wt(u) is enough to capture this concept. In this sense, 



Properties I^J, [5J% are valid for W((u) (with the proportionality factor J^a c <y m the 
second case). 

With regard to the s-e flux vectors, all of them are combined into a single r-parameter 
family of s-e flux vectors which can always be written as 



T t (u) = -TV-.a.^K 1 • • • u^u^ = Y.Ct ( t p? 



T=l 



where the non- negative constants CV are appropriate linear combinations of the c CT . From 
this one sees that Pt(u) is causal and with the same time-orientation than u, so that 



Property 3.7 is also kept. 



In fact, from Property |4.1| , Theorem &A and definition (^g) we have 



Corollary 5.1 The DSEP holds for the general s-e tensor T. 

22 



In other words, a (2r)!/2 r -parameter family of s-e tensors satisfying the good mathematical 
properties mentioned above and including the fundamental DSEP has been constructed. 

A possible way to avoid this freedom (if this is desired) is to take the completely 
symmetric part of the general s-e tensor, which coincides with the fully symmetric part 
of the basic s-e tensor (|): T (MlM .. M2r _ 1Ai2r ) {<} oc T( /il|t2 ..^ 2r _ ljU2r ) {t}. This has no in- 
cidence on the super-energy density, in which definition only the completely symmetric 
part T , (^ 1 ^ 2 .../i 2r _ 1 ^2r) {*} ^ s relevant, and selects a unique s-e flux vector. In fact, in some 
occasions this completely symmetric part arises also in relation with the existence of con- 
served currents constructed from s-e tensors, see Sections |6| and [7]. Nevertheless, I have 
preferred to maintain the generality and the general s-e tensor ( p8[ ) will be used too. 

To illustrate these points, let us consider the simplest cases explicitly. For r = 1 



(simple p-forms XI), the basic tensor (12-13) coincides obviously with the general s-e 



tensor TI\ M \ Eui > = cTx^ I Eu > so that this case is trivial. When r = 2 (double (p, i 

forms ^"[p],[ g ]), the general s-e tensor T a p\^ jiTy u > is constructed from the basic s-e 
tensor fll6Hl7|) as follows 

T a/3 A M {■ K '[p],[ ? ]} = ClT a /3\n { K [p],[q]\ + C 2 T aA/3 M {^[p],[g]} + C 3 T aflX p \ K \p],[q]} + 

+C4T\l3aii \K]p],[q]\ + ^T^\ a ^K^^ + C%T\ mfi |-^[p],[g]| • (29) 

Here ci, . . . ,c% are assumed to be non-negative. This general s-e tensor does not have 
any symmetry in general. However, it is automatically symmetric in the interchange of 
(a/3) <-> (Xfi) whenever the original T a p\^ < -K"[p],[q] \ has this property, in which case the 
following notation will be used: c\ = c\ + cq, &2 = C2 + cq, 03 = 03 + C4, C4 = C5 = cq = 0. 
Thus, there appears a six-parameter (4!/2 2 = 6) family (reduced to a three-parameter one 
when the interchange between pairs holds for T^x^ < ^"[piy \ ) of s-e tensors satisfying the 

fundamental DSEP. Furthermore, T a px^ \ Kw>]M \ ls symmetric in a[3 iff C2 = C4 and C3 = 
C5 (or C2 = C3 in the special case), and symmetric in A/Li if and only if C2 = C3 and C4 = C5 
(respectively C2 = 63). This provides a three-parameter (resp. two-parameter) family of 
s-e tensors with the same symmetries as the basic one. Concerning the interchange of 
pairs, T a pxn \ K\p],[q] \ ls symmetric with respect to the exchange (a/3) «-> (Xfi) if and only 
if c\ = cq and 03 = 04, or in general if the basic one had already this property. 

For triple (p, q, s)-forms -<4[p],[g],[s]> the general s-e tensor ([J8|) becomes a 6!/2 3 = 90- 
parameter family of s-e tensors, and thus the formula analogous to (gjj) is very large and 
will be skipped here. 

6 Application to physical fields. (Super) -energy 
tensors. 

The applications of the above mathematical constructions to the real physical fields is 
considered in detail in this section, containing three different subsections devoted to the 
gravitational field, to typical massless fields (scalar and electromagnetic ones), and to the 
corresponding massive fields (scalar and Proca fields), respectively. Some particular cases 
of interest will be remarked, such as that of General Relativity. The traditional s-e tensors 
of Bel |^, P, Bel- Robinson |4], |?], [S5f| , and Chevreton [24] will be rederived as particular 
cases of the general construction. 
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6.1 The gravitational field 

Let us assume that the gravitational field is described by the curvature tensor R a p\ fl of 
the Lorentzian manifold (V n ,g). Given that R a p,\p is a double symmetric (2,2)-form, the 
basic s-e tensor for the gravitational field is simply the appropriate restriction of (jig): 

TaV \R[2],[2]} = \ (R a pMR^ a + r-^T R * R Pp,W*.*n + 

I ll ' l, J 2 V (n — 3j! ap,\a 4 ...a n 

j f> f J^Pp4---Pn,^o- _, j^ j^/3p4...p n ,fia4...a n \ 

(ra-3)! a P4 :.. Pn M [( n _3)!] 2 a/J4*.pn,Aff4*.<r„ y 

which is a straightforward generalization of the original definition given by Bel [|], y]. 
As remarked at the end of section |3|, after expanding the duals the expression for this 
generalized Bel tensor does not depend on the dimension n, and its form (|17| ) reads (see 



[24j for n = 4, and also J4T 



B a /3\p = Tapx/j, |-R[2],[2]J — Rap,\aR(3 P ' ' p + Rap,pcrR(3 P ' \° ~ 
-7)9oi{iRpT,\aR PT ' p — ~9\pRap,aTRl3 P '' 7T + T.9al39\pRpT,crvR PT,aV (30) 



from where properties ( |18|) are manifest and also 

B a p\p = -B( a/ 3)(A^) = B X pa/3- (31) 

Furthermore, using the first Bianchi identity R a [p\fj,] = if necessary 

B p pX p = B x / P = ^^ (R prM R pT 'p a - \gxpR P r,auR pT ^ , (32) 

(4 - n) 2 

-£> p a = 3 t^pT,avt^ ' j I""] 

o 

1 Ta 1 

B P (3pp = R Pa R/3p,pa — 7;Rf3p,TaRp P ' Ta + o9/3pRpT,auR PT ' aU , (34) 

Z o 

77—4 

£?"V = #*^ + ^Rpr,auR PT ' UU - (35) 

o 

Due to the general results proved in the previous sections, the generalized Bel tensor 
has the DSEP and the generalized Bel s-e density is defined by 

W B (u) = B a ^u a u^u x u^ 

so that 

Wb (u) = \ ([§EE] 2 + [%EH] 2 + [%HE] 2 + [§HH] 2 ) > 

where (§Ee) , (§Eh) , (§HE) and (§Hh) are the electric- 

V" /[1],[1] ^ /[l],[n-3]' ^ A«-3],[l] V" /[n-3],[n-3] 

electric, electric-magnetic, magnetic-electric and magnetic-magnetic parts of the Riemann 
tensor according to the definition of Section 0. For General Relativity (n = 4) they were 



first introduced by Bel in |3[, see also [19]. These E-H parts satisfy in general 



(See) =(*ee) , (§eh) =(§he) 

\ / pv \ /up \ / p,vi...v n ^3 \ / 

(§hh) = (§hh) 



3,pi...p n -3 
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Again, Wb vanishes if and only if the whole Riemann tensor vanishes too 

{3u such that W B (u) = 0} <^=> B a/3 x^ = <==> R a px^ = 0. 

Concerning the generalized Bel s-e flux, from properties (31) and the results of Section || 
it follows that there is only one independent such vector, given by 



pgw 



-B a ^ v vPu p u v 



which is a causal vector with the same time orientation than u. 

In general, from the analysis presented in Section [51 the general s-e tensor for the 
gravitational field is the following 3-parameter family of tensors 



^afiXfi = ClBafiXfj, + C2B a \f3^ + C3B afl Xj3 

whose symmetry properties are 

^>a/3Xfi = ^>/3afiX = ^Xtia/3- 

However, we also have 



(36) 



If c 2 = c 3 

If ci = c 2 = c 3 



a/3Xn — JI5 (a/3)(A/x)) 
®a/3A/i = ®(a/3A/j)- 

and the second one leads to a tensor 



Both cases have been recently analyzed in |9 
considered also in pq| . 

One of the important properties of the generalized Bel tensor is the expression for 
its divergence, which does not vanish in general. By using the second Bianchi identity 
V[ u R aj 3]Xfi = 0, one can deduce from ( pC|) 

\7 a B a ^ 1 = R^J^f + R^ a J Xap - h^R^J^f (37) 

where Jx^p = — J^x/s = ^xR^p — V M i?A/3- Thus, the fundamental result that B is 
divergence-free when the 'current' of matter J\^p vanishes has been proved. More pre- 
cisely 



Theorem 6.1 If Jxu£ = then the generalized Bel tensor ft30[) as well as the 3-parameter 
family of tensors (poj) generated from it are divergence-free. This includes all Einstein 
spaces (i.e. with R^ u = Ag^ u ). ■ 

As is well-known, the Riemann tensor can be always decomposed into the Ricci tensor 
Rp^ = R p j3 Pfl and a traceless part denoted by C a pxfi and called the Weyl tensor |3(j, |53, 5S[. 
This tensor is conformally invariant |36|, 53 1 and possesses the same symmetry properties 
as the Riemann tensor but is also traceless 



a 



a[3,Xfj, 



c. 



a/3], [Am] i 



a, 



IPM 



0, c p 



P>Pl* 







so that C is another double symmetric (2,2)-form. The explicit formula relating these 
curvature tensors is [36, [53, 59] 



Raf3,X/j, — Cap^Xfi + 



n 



Ra[X9ii]/3 — R(3[X9fi]a 



R 



(n — l)(n 



(gaxgpti - ga^gpx) (38) 
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where R stands for the scalar curvature R = R p p . Due to the fact that the gravitational 
field equations usually involve the Ricci tensor in terms of the energy-momentum tensor 
(such as in the Einstein field equations R^ u — (l/2)g flu R = kT^), one thinks of the Ricci 
tensor as the part of the curvature directly related to the matter fields, and of the Weyl 
tensor as the free gravitational field induced by the sources. Of course, one can construct 
the basic s-e tensor for the Weyl tensor, which generalizes the classical Bel-Robinson 
tensor ||, 78, ^, 97] constructed in General Relativity (n = 4): 



T a V {C[2],[2]} = \ (c apM C^° + -^— C . C tol"*™ 9 » + 



[jl — o)\ ap,\tJ4...<j n 

1 * 1 * * \ 

(n-3)! a Pi :.. Pn M [( n _3)!] 2 a p£.. Pn M*-an j 

This can be re- written, analogously to (|30|), independently of the dimension n as (compare 



with [29 



%x/3X/j. = T a 0\p (C[2],[2]| — C aP) XaC/3 P, p Cr + C ap ,p,oCp P "\ — 
— -^9apC pT ^\ a C pT, p CT — -g\ p C ap ^ T Cf3 P,CJT + -g a (39\pCpT,aiyC pT ' aU (39) 

form where one deduces the symmetry properties 

%x/3Xp = 2(Q,/3)(A/i) = TxpaP- (40) 

Taking into acount that the Weyl tensor is traceless, the formulae analogous to (p2[]35|) 



arc 



* pXp — *\p, p — ^ I ^pr,Xa^ 'ju ~~ ~:gXp^- J pr,av^ ' I , (41) 

(4 - n) 2 
TP P a ° = B C pT , au C pT ^, (42) 

o 

T p (ipp, = -- ( C / 3 PiTCT C /i p ' rCT - -g(3pC pT ^ u C pT ' au \ , (43) 

r? — 4 
TP V = -^Cpr l( r,C^ w . (44) 



The generalized Bel-Robinson tensor ( |39| ) satisfies the DSEP, and one can define the Bel- 
Robinson s-e density Wq- (u), the flux vector Pj- (u) and the corresponding E-H parts of 
the Weyl tensor. All these objects have the same good properties as those considered in 
the case of the generalized Bel tensor, and their definitions can be obtained by simply 
substituting the Weyl tensor for the Riemann tensor everywhere (and T for B). Similarly, 
there is a 3-parameter family of s-e tensors for the Weyl tensor with properties analogous 
to those of B. Nevertheless, all these general properties degenerate in the case of n = 4 
(including General Relativity), where the Bel- Robinson tensor has special properties to 
be shown presently. 

Concerning the divergence of (|39"1), by using (J3^) , the formula for V\ v C a mxp. deduced 
from it and 

n — 3 / 1 \ 

one arrives to 
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Proposition 6.1 If R^ u = Ag^ u , then the generalized Bel-Robinson tensor (31) as well 
as the 3-parameter family of tensors generated from it are divergence- free. ■ 



Obviously, the generalized Bel-Robinsor tensor (|39| ) is always a part of the generalized 
Bel tensor (|30|), as follows from pgj). Actually, by introducing ([38]) into (^) one gets the 
canonical decomposition of B as 

B a f3\ii = XypXfi + MafiX^ + Q a 0Xfi 

where M a pxp is called the pure matter gravitational s-e tensor (it depends exclusively on 
the Ricci tensor) and is given by 

M af 3Xp = T a/3 x^ [R[2],[2] ~ C[2],[2] j , M a 0X/J. = -M(a/3)(A M ) = Mxfj.a/3 (45) 

and Q a f3X[i is the matter-gravity coupling tensor, containing all the terms with products 
of the Weyl and Ricci components and whose explicit expression reads 



Qa0Xfi = ~, — 3^7 ^C 7 (X^)( a R/3)a + 4C" 7 ( Q/ g)( Aj R At)(T + 2R pu [C ' a p (A ° g y,) p - C X P // g a + 

R 

+C/3 P (x Cr 9fi)a - Ca P p a gXn) , Ra/3 = Rf3a = Ra/3 ~ ^7 — ~jTSa/3 

from where it is easily derived that 

2(4 — n) 

Qa/3Xn = Q(a/3)(Xfi) = Qx^a/3, Qa(0X/i) = °) Q P pXp = ~, 7^^° CpXvfJ,- (46) 



n 



The form of Q is in fact simpler in the case of n = 4, where 6 Q a pxp = R(C a x0p + C afJi px), 
as can be proved using spinors HlOOH or traditional methods [pLSf l . Using (|46|) , in general n 
we have 

B a(0Xfj.) = %x(J3Xn) + MatfXfj,) (47) 

so that the completely symmetric part of the generalized Bel tensor decomposes simply 
as the sum of the completely symmetric parts of the generalized Bel- Robinson tensor plus 
that of the pure-matter gravitational s-e tensor. 

The pure-matter gravitational s-e tensor ( |45| ) is interesting because only the Ricci 
tensor of the spacetime is involved, which in most cases is equivalent to saying, due to 
typical field equations, that only the energy-matter content is involved. Given that Ai is 
the basic s-e tensor (||) for the double symmetrical (2,2)-form (i?[2],[2] — C[2],[2]h it follows 
from the results of Sections y and || that the pure-matter gravitational s-e density Wm {&) 
is positive definite, that the pure-matter gravitational s-e flux vector Pj^x (u) is causal and 



that in general Ai a px^ satisfies the DSEP. Moreover, one has from Property |3,4| and 



{3u such that W M (u) = 0} ^=> M Q% = <^=4> R a p,Xv. = C Q/3jA ^ <==^> R^ = 0. (48) 

Even though these properties of M. indicate that it is a good s-e tensor for the Ricci 
tensorjj, one should bear in mind that this is not a s-e tensor for the particular field 
(or type of matter) which generates the gravitational field through the field equations. 



4 Notice that another similar basic s-e tensor can be constructed using that the Ricci tensor R[i][i] is a double 
symmetrical (l,l)-form, i.e., T a fjx^ {R\i]\i\ }- This tensor has all the above good properties too, and in fact it 
can be related to M. a px^i see for instance p9| . 
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A good s-e tensor for a matter field must have a well-defined sense in the absence of 
gravitation (for instance in Special Relativity), and this is not true for A4. Somehow, 
the term gravitational pure-matter s-e tensor has been systematically used to remark that 
this is a gravitational super- energy, even though it appears due to the matter contents of 
the spacetime exclusively. The good s-e tensors for physical fields will be constructed in 
the next subsections. 

From the second property shown in (|46|), or directly from ([47|), it follows that the 
matter-gravity coupling tensor Q does not contribute to the generalized Bel s-e flux vector. 
More precisely 

Proposition 6.2 The generalized Bel s-e flux vector is the direct sum of the generalized 
Bel-Robinson s-e flux vector and the pure-matter gravitational s-e flux vector. 



Proof: Contracting ( [471) with the unit timelike u thrice one gets 

Pb {u) = Pt (u) + P M iu) 
which proves the result. ■ 

Corollary 6.1 The generalized Bel s-e density decomposes as the simple sum of the gen- 
eralized Bel-Robinson s-e density plus the pure-matter gravitational s-e density: 

W B (u) = W T (u) + W M (u) . (49) 



Corollary 6.2 In general n, one has 

Proof: If R^y = the result is trivial. Conversely, if B a p\^ = T a pxn, then Wb (u) = 
Wt (u) so that from (f49[) we arrive at Wm ($) = 0- Using (|48|) the corollary follows. ■ 

Therefore, in Ricci-flat spacetimes there is no need to distinguish between the gener- 
alized Bel and Bel-Robinson tensors. 

The properties of the different s-e densities allow to compare the relative strength of 
the Riemman, Weyl and Ricci tensors in a given (V n ,g). Thus, for instance, one can define 
the three positive scalars 

_ W M (n) _ W T (u) __ W r (u) 

W B {u) W m (u) W b (u) 

relative to any observer u, which are not independent in general and such that, due to 

< qi < 1, < q 2 < oo, < g 3 < 1. 

These may serve to measure the strength of the Ricci tensor in front of the pure Weyl 
tensor or of the whole Riemann tensor. Thus for example, qi can be used to check the 
so-called Penrose's Weyl tensor hypothesis |79| , pH] and similar conjectures concerning the 
entropy of the gravitational field pOJ ]. On the other hand, q\ may be a 'quality factor' 
for approximate solutions of some field equations as explained in [||, 19|. In this sense, 



2S 



notice that eft = if and only if R^ = 0, which is usually equivalent to the absence 
of matter-energy content. Thus, the closest the adimensional number q\ is to zero for 
any given metric and some u, the best resemblance it has to an exact 'vacuum' solution 
for that observer u. Observe that one does not need to know the explicit exact solution 
one wishes to approximate in order to compute q%, so that the goodness of approximate 
solutions may be measured without ever knowing the form of the exact solution. Finally, 
53 gives the relative strength of the Weyl and the Riemann tensors or, in other words, 
of the free gravitational field with respect to the total one. The number q% vanishes if 
and only if the metric is conformally flat, so that q^ measures the departure from this 
condition somehow. This can have relevance in, for instance, General Relativity and other 
typical theories, where the standard cosmological models are conformally flat. 

In General Relativity (n = 4), the Bel-Robinson tensor has been widely studied and 
used, and it is well-known that it is a completely symmetric and traceless tensor ||, 0, [lj, 
1|, |I|, H§, 0, |67|, ||, [74], [78], pa ||, H|, ||, |97|, p], [TOOR. These properties do not hold in 



general, and in fact (see also |4l| ): 

Proposition 6.3 The generalized Bel-Robinson tensor is completely symmetric if and 
only if n = 4, 5. 

Proof: From (]42|) and (H) one gets 

TP a _ TP a — ±2 ~ n )\^ ~ n > i~i ptpT,ov 

* p a ■* pa — ^pr,au^ 

o 

and this must be zero if the tensor T is completely symmetric. Thus, a necessary condition 
for the complete symmetry of the generalized Bel-Robinson tensor is that n be either 4 
or 5. To see that it is also sufficient, let us consider both cases separately. 
Case n = 5 . From formula ( |A.12| ) in the Appendix, if n = 5 

c * x ™cP pa . = AC ap Xa C p P a - 

apa pvrv r r 

—2o^C a p TU C — 2o a C (-^ppru + o a d^Cp Ttau C ' (5U) 



but from the identity (Eh or directly from (A. 13) we also have 



c . Xtu = -9d Ar fQ „<r 

apa y a P °1 

so that using repeatedly the trace-free property of the Weyl tensor, a straightforward 
calcualtion leads to 

* * 

(-1 Xrv^Ppa _ Afi Xa^Pp i q(~i fipf-iXo , a(~i pas-tXfi , 9 c/3 ^Xpru s~i , 

aoa u*u ~ ~ a P pa ~r S^aa l> pp ~T 'i^ap ^ pa ~r z «a ( - / ^ppru ~T 

4-0 X^C rtPpTV _ AXfir 1 n~Xprv _ ArX(~itipTV(~i A-X^xPn (~tpT,av (r-i\ 

-|- ZO ^U a p T v^-< ^O^^apru^ *°a L/ ^ pprv *C < a O ^O pr ,<n/0 ■ K 01 ) 

As ( |50| ) and (^) must be equal, the following identity follows (see also J65[ ) 

Tf „ _ tc v 9 /'/^ Xa/~if3p f-i f3a^<Xp \ (-< pa (~iXji r\(~t f~i/3pTU 

11 IL — O f Z I ^ap ^ pa ^ap ^ £«7 J ^ap ^ per ° p K - J a.prv K -' 

-biC^ TV C WTV + ^CV^™ + ^"C^ + i ($£<$ - <£$) Crr^C"*" = 0. 
But due to Cq,^] = 0, 

Cp&sy\{3 o/^f ( s-y(3<j\p /~i/3p\a 

ap ^ per — ^^appa I ^ ^ 
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so that the previous identity, using ( 3£ ) , can be simply written as 

if n = 5, 



1~ag\[i — %t\P/x — 



which together with 
Case n = 4 (see j 



40) proves that T a pxp is completely symmetric in n = 5. 
. From ( |A.12| ) in the Appendix, if n = 4 



CAfjp/3p /^f Act /~ij3p X^r* 1 /~ij3prv rft /-~i\pTV s~i , xft X^f* 1 r~ipT,av 



op 



pO" 



M°" „^ v - / apTi/i-' r ) u a , - y ^pprv T * u a u p^ pr.av^ 



but the Lanczos identity ( A. 15 ) implies 



C 



op 



-c 



Ar 



op 



(52) 



so that we obtain 






J^aprv 



fipru 



J cc J \x KJ PT.pv^ 



or equivalently 



C a/)ri ,C^ - ifc,C^ = if n = 4. 



J a^pr 



Observe that this implies the complete tracelessness of the Bel- Robinson tensor in n = 4 
due to flU). From (||) it follows that 



C * = C . =C a/3A p if n = 4 

op,Ap op,Ap 



(53) 



so that only one independent dual can be constructed in ra = 4. Thus, the traditional 
formulas for the Bel-Robinson tensor in n = 4 can now be written [Q, 0. |78|, Q^, 97, p9| : 



T Q V = C apXu C^+ Capx.6 0PIMT , if n = 4, 

^/sAp = C ap \ a Cp p '[/ + CappcrCi3 p \ e7 - -g a f3gxpC pT! auC pT,cr ' / , if n = 4. 

o 
Using C a r^^i = one can write 



(54) 



L'apAtr — (,^«pAcr ~r L^ApcwT ~r ^aXpaj 



and analogously 



CvapAcr „ ^ CvapAcr i C/Apacr "T U ctXpa 



so that formula d54f) becomes 



^a/3Ap — T {CapXa + Ca P oo-) {C^ P ^ + Cyfp ) + -C a \ pa Cj3p PCr + 



(sapXa t C/Apocr I I O/J p t C/p /3 1 . C/aApcrC'p'p 



1 

4 



if n = 4 
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and using here that, from (A. 7) or directly from (Q), 

C aXp(T C^ pu + CaXpaC^ pa = if n = 4 

one obtains an expression for T a p\^ which has manifestly the property 



%x/3Xn 



T> 



\f3ap, 



if n 



which together with ( |40| ) proves the result. ■ 

From the above proof, only one independent dual of the Weyl tensor is possible in 
n = 4, see (|53|). Thus, the E-H decomposition of the Weyl tensor in n = 4 is very simple, 
because from (52) and (|53"|) it follows 



f $EE) 



C uHH) 



C 



EH 



\n> 



$he) 



if n 



fJ,V 



which are called the electric and magnetic parts of the Weyl tensor relative to u if n = 4 
HI H> 0i H3' 15' —\- These tensors are spatial relative to u, symmetric and traceless; 
they were first introduced by Matte [70]. An analysis of the positivity properties derived 
from the DSEP for | 



appears in |12[] . Such an analysis, and its consequences, may be 
translated to the general (g) in principle. 

A traditional question concerning the Bel-Robinson tensor is what type of physical 
quantity, if any, is described by it. In this sense, one must notice that the physical dimen- 
sions of T or B are L~ 4 , where L means length. This has been usually separated into two 
L~ 2 's, and led to two different classical interpretations: either both L~ 2, s are on the same 
footing describing an energy density each so that Wb (u) is an energy density "square" 
(]T3|, |1S| , |34j and references therein), or one of the L~ 2, s is a 'true' L~ 2 and Wb (u) repre- 
sents energy density per unit surface (or second derivatives of energy density). The first 
interpretation was supported because, for instance, in Einstein-Maxwell spacetimes and 
n = 4, the pure-matter gravitational tensor takes the simple form M a p\^ = RapRx^, 
which through Einstein's field equations can be rewritten as the square of the energy- 



momentum tensor [19]. This led to the possibility of finding a 'square root' of the Bel- 
Robinson tensor, which would be a good two-index tensor with appropriate physical di- 



mensions. Unfortunately, this square root exists only in some particular cases [19, [l 



A generalization of the square root to two different factors is always possible [19, 13], 



but in this case the factors are not unique in general [19, 13]. Thus, the first interpre- 
tation seems dubious. On the other hand, the second interpretation seems correct as 
can be deduced from several independent results p2| , 

i Hi, n, m, mm n, 
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Notably, the relationship 



57, B6 9 J 



E a = const. W r (u) \ r=0 a 5 + 0(a 6 



between any of the quasilocal energies E a of small 2-spheres of radius a in vacuum and the 
Bel-Robinson s-e density Wq- (u) supports clearly this view (u orthogonal to the 2-sphere 



and n = 4). Similarly, the results in [42] point into this direction and, more importantly, 
Teyssandier |95| has recently found a formula for the 'super-Hamiltonian' associated to 
the super-energy of the scalar field (see next subsections) in Special Relativity from where 
one checks that this super-energy is interchanged in quanta of fnu^/c 2 , where u>k is the 
frequency of the fc-mode. All in all, due to the independence of these three results, it 
seems that the second interpretation is valid, at least in n = 4, as also suggested in J3(| . 
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To end this section, let us finally make some comments about the higher order super- 
energy tensors for the gravitational field. Instead of considering the Riemann tensor as 
fundamental field, one can also take its covariant derivative and repeat all the construction 
above. Specifically, one can consider the tensor V u R a f3Xp as a triple (l,2,2)-form and 
construct its corresponding s-e tensor (||), T a p\^ TU \ Vfiiiipip] f > an d analogously for the 

Weyl tensor, etcetera. This provides an infinite set of (super) -energy tensors for the 
gravitational field, one for each natural number k > 1, according to the following 

Definition 6.1 The ( super ) k -energy tensor of the gravitational field (k > 1) is the tensor 
(j^j for the (k + l)-fold (1, . . . , 1, 2, 2) -form describing the (k — l)-th covariant derivative 
of the Riemann tensor, that is 

T\ lfll ...\ k n k \ k+1 n h+1 < V[i] . . .V^ -R[ 2 ],[2] 

I fe-l 

A similar definition can be given using the Weyl tensor. Prom Ql9D, the explicit formula 
for the (super) 2 -energy tensor of the gravitational field is 

T a /3Xpru |V[l]i?[2],[2] j = V Q ,i?ApT(rV / gi? At P !/ Cr + V pRxprv^ aRp P v" + V ' a Rpp,Ta^ ' pR\ f \° + 

+V a R Xpua Vf 3 R/ T CT - g af3 (v^Rxpra^R/^ + V c i2 Ap ^V^ 

— ydXfJ, [y aR^fyra^ pR u a + V ' a R^p Ua S/ ' pR 'f P T °' 

- X -9tv (v a Rx P <;*VpR/^ + S/ a Rp p(a VpR X P ^) + 

+\g a pgx f yc R iprMR^u' T + ^gapgruVtRxp^R/^ + 
+\gxpgruS/ a R (py(7 V p R^ a - papgxpdrvVsRcrreV 6 &"'-(&) 

These (super) -energy tensors may have some relevance at points of V n where the Riemann 
tensor vanishes but such that some derivative of Riemann is non-zero there (so that every 
neighbourhood of the point has gravitational field). From the general results of sections 
H] and U it follows that all the (super) fc -energy tensors satisfy the corresponding dominant 
(super) fc -energy property, and one can define the (super) ^-energy densities relative to u. 
Moreover, we have 

Proposition 6.4 The gravitational ( super ) k -energy (tensor) vanishes in a domain D C 
V n if and only if the (k — l)-th covariant derivative of the Riemann tensor is zero in D. 
In particular, all gravitational ( super ) k -energy tensors vanish in flat regions of (V n ,g). 
If the gravitational (super) k -energy (tensor) vanishes in D, so do all the ( super ) k -energy 
tensors with k > k. The (V n ,g) of constant curvature are characterized by the vanishing 
of the gravitational (super) 2 -energy (tensor). ■ 

Thus, for example, in Special Relativity the gravitational (super) -energies (and therefore 
the corresponding (super) ^-energy tensors) vanish. Besides, the previous Proposition 
clarifies somehow the reasons why the Bel-Robinson and Bel tensors arise naturally in 
General Relativity: from the equivalence principle, the energy density of the gravitational 
field can be always set to zero at any point x € V n by appropriate choice of the observer, 
and thus the natural concept at the given point is the super-energy density. Further 
related considerations regarding these tensors will be made in the next subsections. 
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6.2 Massless fields 

General massless fields can be considered from the point of view of super-energy tensors 
and, sometimes, the use of spinors simplifies things substantially [12, 13 , 14, [0], 24, ^4, 8C[. 



However, the spinor simplification is usually restricted to dimension n = 4 and thus, in 
this subsection, only the outstanding cases of the electromagnetic and the massless scalar 
fields are going to be treated in general dimension n. 

6.2.1 The massless scalar field 

Consider a massless scalar field 4>, that is, any function </> satisfying the n-dimensional 
massless Klein-Gordon equation 

V p W = 0. (56) 

This is usually referred to as a minimally coupled massless scalar field, and is not con- 
formally invariant in general [p0| . As the basic object is </>, the simplest s-e tensor (H) 
for (f) would simply be T{(p} = 4> 2 /2. This object is obviously positive but it seems to 
have no physical relevance. However, one can also take the first covariant derivative V»</> 
as starting object so that a case of physical interest arises by constructing the tensor (g) 
associated to Vhi0. Using (|l3|) one gets 

Tx^ [1]( f>} = Va<£V^ - IgxpVpWh V{V[i]0} = (l - |) V p 0W (57) 

which is, in fact, the standard energy-momentum tensor of a massless scalar field. The 
tensor (E^) is symmetric and identically divergence-free if the field equation ( |56| ) holds. 
These two properties lead to the existence of conserved currents for the scalar field in any 
spacetime with a Killing vector field £, that is, with a vector field satisfying [36, ^] 



V M &, + V„£„ = 0. (58) 

The well-known idea is to construct the current 

which due to fl58| ) and the properties of Tx M {Vn]0} satisfies 

V^ (V^^l) = T^IV^^V^) = 0. (59) 

Notice that in this simple calculation both the divergence-free as well as the symmetry 
properties of T p/i {Vri]<^} are needed. Once that we have a divergence-free current such 

as j M (Vmi(^;£j, conserved quantities are readily constructed by means of the Gauss the- 
orem |53| and the integration over appropriate domains of the manifold V n (sometimes 
called Gaussian flux integrals, 0]). In flat (V n ,g), they provide n(n + l)/2 independent 
conserved quantities for the scalar field (one for each independent Killing vector) related 
to the energy, the momentum and the angular momentum of the scalar field. In general 
non-flat (V n ,g), they give some conserved quantities related to the intrinsic symmetries 
of the spacetime, and to the corresponding combinations of the pertinent quantity for the 
scalar field together with the gravitational field (described by the curvature of the (V n ,g)). 
Observe finally that only a conformal Killing vector £, that is 
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is necessary in case that the tensor T\^{V[i]4>} be traceless, because then j^ ( Vm^C) is 
also divergence-free in general. In the case of the scalar field, this only happens for n = 2 



as can be seen from (57) 



One can go on with the above construction and use now the double symmetric (1,1)- 
form V a V/s4> as the starting object, so that the corresponding tensor (||) becomes by using 

(0) 

S a p\n = r a/3A M {V[i] V[i]0} = V a V \4>V ^V ' p(j) + V Q V ^V xV p(j) - 



In fact this tensor was previously found by Bel |J and Teyssandier [95] in Special Relativity 
(n = 4). In general, it will be called the basic s-e tensor of the massless scalar field. Its 
symmetry properties are 

S a 0\n = 5'( Q/ g)( AM ) = Sxfj.a/3 

so that, following the considerations of section [|, the general s-e tensor for the massless 
scalar field is a three-parameter family given by 

^apXfi = biSapXfj, + &2<SaA/?/i + bsS ap xp 

for arbitrary non-negative constants b\ , 62 , ^3 • The symmetry properties of § are the same 
as that of B 

^apXfi = ^f3afiX = ^Xfiaf3 
If 62 = 63 => Sa/3A/i = §(a/3)(A M ), 

If 61 = 6 2 = 63 => §a/3A/i = ^(a/3A^) • 



The last two cases have been also found in 95] for the case of Special Relativity (n = 4). 



All in all, the tensors S a pXf_i and E> a pxfi have mathematical properties analogous to 
those of the generalized Bel tensor and the general s-e tensor of the gravitational field. 
This analogy can also be confirmed on physical grounds by several independent methods, 
such as a) the analysis performed in |95| ], b) the interplay between the propagation of 
discontinuities of the Riemann tensor and that of the second derivatives of <p leading to 
mixed conserved quantities along null hypersurfaces |]3S[ ], or c) the existence of general 
conserved currents involving necessarily both the generalized Bel tensor and the tensor 
(|60|), sec [90 ] and section |7[ Therefore, it seems that the tensor S a px p represents for the 



scalar field the same type of physical quantities as the traditional Bel and Bel-Robinson 
tensors do for the gravitational field. Thus, the super-energy density and the s-e flux 
vector of the massless scalar field can be defined using the tensor ( p0|) . 

Concerning the divergence of S a pxu, a straightforward computation using the Ricci 
identity J36J, ||, 0] leads to 



-V^ (2V<V (A ^ )p/ 3 + <?A M ^ T V"V>) (61) 

so that, as we can see, all the terms on the righthand side involve components of the Rie- 
mann tensor. Notice that, due to the symmetry properties of S a px p , only one independent 
divergence can be computed. Thus, one has the important result: 
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Theorem 6.2 In any flat region of(V n ,g) (that is, with vanishing Riemann tensor), the 
super- energy tensor of the massless scalar field f p^ ) as well as the general 3-parameter 
family of tensors generated from it are divergence-free. ■ 

In other words, in absence of gravitational field (represented by the Riemann tensor), the 
general super-energy tensor of the scalar field is divergence-free. In fact, this leads to 
conserved currents and quantities for the scalar field in flat (V n ,g) in the same way as was 
shown before for the energy-momentum tensor of <j), given that any flat (V n ,g) always has 
n(n + l)/2 independent Killing vectors |j(]]. Thus, define 

for any given Killing vector £, where the symmetry properties of ( pOJ ) have been used to 
write the last expression. This j is divergence-free in flat (V n ,g) because then 

V M ^ (V[i] V M 0; |) = S (w<Tr) V M (£ p Ut) = in flat (V n , g). 

Actually, one can put different Killing vectors in the definition of j if at the same time 
only the completely symmetric part of (pCf) is used. Thus, the vector fields 



U (V[i] V [!](/>; II, 6, 6) = S^^tftftf (62) 

are divergence- free in flat manifolds for arbitrary Killing vectors £1 ,£2 ,&■ This provides 
a total of 

( N + 2 \ 

[n-i) 

conserved currents (|62| ) in flat spacetimes, where iV stands for the number of independent 
Killing vectors: n{n + l)/2. These are currents quadratic in the second derivatives of 
(j). Then, one can find the corresponding conserved quantities using the Gauss theorem, 
as always. In the particular case of n = 4 there are 220 currents in the set (^), to be 
compared with the results found by Teyssandier |95| ], 

Continuing the process, one builds the (super) 2 -energy tensor T' Q ,^> iMT ^{Vr 1 iV[ 1 ]V[ 1 ]0} 
associated with the triple (l,l,l)-form VqV^V^, which using ( |l9| ) is given by 



Tafix^A^ [l^ll^li]^} = V Q V A V T 0V /3 V M V^ + V /3 V A V T 0V a V M V ! ,0 + 

+V Q V M V r (/> V^V A V^ + V Q V A V v V /3 V M V r - 

-9*0 (V P V A V T V%V^ + V P V A V^ V V M V T 0) - 

-0A,. (V«V p V t V/jVVv^ + V a V p V^V^V^V^) - 

-5™ (V Q V A V p V/jV^V^ + V a V^Vp0 V/jVaV^) + 

+5a/35A M VpV CT V T VVV,^ + 9 a /39r^V p V A V CT V%V> + 

+9A M 5r, V Q V p V CT V^ VV> - -g a 0gx„9rv V 7 V p V CT </> V^VV> (63) 

and so on. This produces an infinite set of basic (super) fc -energy tensors, one for each 
natural number fc. 
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Definition 6.2 The (super) -energy tensor (k > 0) of the massless scalar field is the 
tensor ftQ) for the (k + l)-fold (1, . . . , l)-form describing the (k + l)-th covariant derivative 
of 4>, that is 



?A 1M1 ...A feWb Afc + i«b+i { V [i] • •• V [i]^ 

fc+i 

Each of these tensors describes the same physical characteristics of the massless scalar 
field as the corresponding gravitational (super) ^-energy tensor does for the gravitational 
field. They may again have some relevance at points of V n where the derivatives of a 
given order, but not of all orders, of the scalar field vanish. From the general results of 
sections and B the dominant (super) fc -energy properties hold, and one can define the 
(super) -energy densities relative to u. Moreover, the following fundamental result holds: 

Proposition 6.5 The ( super ) k -energy (tensor) of the massless scalar field vanishes at 
a point x € V n if and only if the (k + l)-th covariant derivative of cj) is zero at x. In 
particular, all ( super ) k -energy tensors of the massless scalar field vanish in a domain 
D C V n if (j) is constant in D. If the (super) -energy (tensor) of cj) vanishes in D, so do 
all the (super) -energy tensors with k > k. The manifolds containing a constant V(j) are 
characterized by the vanishing of the super-energy tensor of (f). ■ 

Notice that the (super) -energy tensors of the massless scalar field have a precise and 
definite meaning in flat manifolds (V n ,g), and one wonders whether these tensors give also 
rise to conserved currents similar to (p2|). In general (V n ,g), one can compute the diver- 
gence of the (super) fc -energy tensors to obtain very long formulae involving the Riemann 
tensor. Observe that now there may appear, in fact, several different divergences for these 
tensors, depending on the index used to contract with the covariant derivative. Neverthe- 
less, in flat (V n ,g) it is obvious that Vm • . . Vm</> is completely symmetric, and thus all 
the (AYAir)-pai r s of the (super) -energy tensors can be interchanged in this case. In other 
words, in the case of flat (V n ,g) only one independent divergence of the (super) ^-energy 
tensors exists. Computing it in the special case of flat manifolds one straightforwardly 
gets using (pq) 



VAir Ai W ...A*M*A*+iA**+i ,j V[i] . . . V[i] (f>\ = in flat (V n , g). 

fe+i 



Theorem 6.3 In any flat region of (V n ,g), all the basic ( super ) k -energy tensors of the 
massless scalar field, as well as the general (2k +1)!! -parameter family of tensors generated 
from it according to pjj), are divergence- free. ■ 



( ?2 • • • ?2fc+l?2fc+2 



From this theorem one is led to construct the following currents 

jx | V [i]- - -V[i]^ ;6. • • ■ >6fe+2 = T { x fH ...\ h+1 ix h+1 ) \ V[i] • • • V[i] # < 
fe+i / I fe+i 
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which are divergence- free in flat manifolds (V n ,g) for arbitrary Killing vectors ^2> • • • > £,2k+2- 
This provides an infinite set of conserved quantities in flat (V n ,g) for the massless scalar 
field. This result is not surprising, because if (f> satisfies the Klein-Gordon equation ( |5(j| ) 
in flat (V n , g), then the scalars V| ai | . . . V| a k+1 \4> f° r fixed values of a±, . . . , a^+i so do, 
and their respective "energy-momentum" tensors give rise to an infinite set of conserved 
quantities. For n = 4 a clarifying discussion can be found in [p5|. 



6.2.2 The electromagnetic field 

Let F pv = Ft pu i be an electromagnetic field satisfying the source-free Maxwell equations 

F„ u = V P A V - V V A P ( => V [T F H = ), VpF"" = (64) 

where ^4 M is the electromagnetic potential, which is defined up to the gauge transformation 
A p — ► A p + V^y? for any arbitrary function ip. One could first consider the electromag- 
netic potential A as starting tensor, and build up the tensor (^) associated to Am , given 
by 

T v {A (1] } = A x A p - \gx»A p AP. (65) 

However, this tensor is not invariant under gauge transformations and thus is not well- 
defined. Nevertheless, this tensor will have some role to play in the case of the massive 
Proca field, to be studied in the next subsection. 

One thus must take F, which is gauge-invariant, as the tensor describing the electro- 
magnetic field. Constructing the expression (J6|) associated to F, or using directly fli~3|), 
one gets 

T Xp {F [2] } = F Xp F/ - \g^F pa F^, T/{F [2] } = (l - ^ F p(T F^ (66) 

which is the standard energy-momentum tensor of the electromagnetic field. This tensor 
is symmetric and identically divergence-free if the Maxwell equations ([34|) hold, but it is 
traceless only in n = 4, as seen from (|66|) . As in the case of the massless scalar field, this 
leads to divergence- free currents whenever there is a Killing vector £, defined by 

f (F [2] ; |) = i p T^ {F [2] } => V,f (F [2] ; £) = 0. 

In (V n ,g) of constant curvature, there are n(n + 1)/2 independent such currents, which are 
related to the fluxes of energy, momentum and angular momentum of the electromagnetic 
field. In general non-flat (V n ,g), they provide conserved quantities arising because of 
the intrinsic symmetries of the spacetime, and related to combinations of the pertinent 
quantity for the electromagnetic field together with the gravitational one. Notice also that 
a conformal Killing vector £ will give rise to a divergence-free current j^ [F[ 2 y, M in n = 4, 
when the energy-momentum tensor is traceless. Thus, in conformally fiat V4 manifolds 
there appear 15 such conserved currents. 

In order to get the basic super-energy tensor of the electromagnetic field, the double 
(l,2)-form Vq-F^i, can be used as the starting object, so that the tensor @ becomes by 
using (|T 



E a px p = T a p Xfl {V {1] F [2] } = VaFxpVpF/ + V a F^VpF x P - g a pV «F Xp V° F/ - 

-\gx^ a F ap V^P + ^g a pg x ^ T F ap V T F^ . (67) 
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Its symmetry properties are 



E a /3\n — E (a0)(Xfi) 



but it is not symmetric in the exchange of a/3 with A/i. Thus, following the considerations 
of section ||, the general s-e tensor for the electromagnetic field is the six-parameter family 
(g|) given by 

^ctPXfj, = ClEapX/j, + C2E a xi3fj, + C 3 E afl x/3 + Q^A^ + C^E^\ a + CqEx^P (68) 

for arbitrary non-negative constants c\ , C2 , C3 , C4 , C5 , cq . The tensor E has no index sym- 
metries in general. However, there are the following particular cases: 

If C 2 = C4, C 3 = C 5 => E a/3A ^ = E( Q/3 ) Am , 

If c 2 = c 3 , c 4 = c 5 => E a/ aA/i = E a/3 ( AAt ), 

If Cl = C6, C3 = C4 => E Q/3Am = E AmQ/3 , 

If c 2 = c 3 = c 4 = c 5 =>■ E Q/3Am = E( a/3 )( AM ), 

If C 2 = C3 = C4 = C5, Cl = C6 =>■ E Q/3Am = E( a/3 )( AM ) = E AmQ/3 , 

If Cl = C 2 = C 3 = C 4 = C 5 = C 6 => E q/3Am = E( a/ 3 A/i) . 

The last but one case was considered many years ago by Chevreton [24] in n = 4, see also 



], and this together with the last case have been recently considered in |95[] in Special 
Relativity. 

The same considerations as those made for the case of the massless scalar field apply 
now, and thus the tensors (^) or ( |68| ) seem to describe the same physical quantities 
for the electromagnetic field as the traditional Bel and Bel-Robinson tensors do for the 
gravitational field, or the tensor ( |60| ) does for the massless scalar field. Thus, the super- 
energy density and the s-e flux vectors of the electromagnetic field can be defined using 
the tensors (|67|) or (|6§|). 

Due to the symmetry properties of E a px/j,, one can compute two independent diver- 
gences, which after a straightforward calculation using (pih and the Ricci identity become 



VpE Pp\n = 2 (R p /°F T ^ + RP aT ^F r p\ V CT F A ^ + 
+2V /3 F<% (R x ^F a p + F x \R ap - R x ^ av F av ) + 



+g x » [R^ T F up V T F^ - R Ta F T P\7' 3 F% - -R^F^F^ 
and 



v E a/3 P » = 2 ( R °(<* F _ ! F °t R (* \ v /3) F w + 3y(<*F rP)t[^ F p] t + 



_ g afS 



R™F ap - l -F° T R v pa ^\ V v F»p + ~V v F op Rrb"Frt T 



As is explicit, these righthand sides vanish if the Riemann tensor is zero. 

Theorem 6.4 In any flat region of(V n ,g), the basic super-energy tensor of the electro- 



magnetic field mFa ) as well as the general 6-parameter family (68) generated from it are 
divergence-free. ■ 
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Thus, in absence of curvature, the general super-energy tensor of the electromagnetic field 
is divergence-free. This again leads to conserved currents for the electromagnetic field in 
flat (V n ,g) similarly to those shown before for the scalar field. Thus, the vector fields 



Jn (V[i].F[2];£l,£2>£3J = E^ p(TT ^i p ^2 a S,3 T 

are divergence- free in flat manifolds for arbitrary Killing vectors £i,£2>£3- The total 
number of independent such conserved currents is the same as that given for the scalar 
field. Again, these are currents quadratic in the first derivatives of F and one can find 
the corresponding conserved quantities using the Gauss theorem. 

As in previous cases, higher order s-e tensors can be constructed. The basic (supers- 
energy tensor of the electromagnetic field is the expression (||) associated to the triple 
(l,l,2)-form V p V a F pu . Its explicit expression is 

T a ^ TV {VflVflFp]} = V Q V A Fr P V^V M F/ + VpVxFrpVaV^FS + 

+V a V M F Tp V /3 V A F I / + V a V x F up V p VpF T P - 
-g a p (V VxF T/ >V°VpF v t> + v a y x F vp v a v^F T p) - 
-gx» {V a V a F Tp V p V a F v P + V a V 9 F vp VpV a FS) - 

-\gru (V a VxF ap V/3V p F^ + V a V p F ap Vf3X7 x F°P) + 
+5 a /35A M V CT V 7 F rp V CT V 7 F/ + ^ Q/ 35r,V 7 V A F CTp V^V M F^ + 
+^A^r,V a V 7 F ap V /3 V^F CT " - ^gapgx^r^sVyF^VF^. (69) 

In general we have 

Definition 6.3 The (super) -energy tensor (k > 0) of the electromagnetic field is the 
tensor ^) for the (k + l)-fold (1, . . . , 1, 2)-form describing the k-th covariant derivative of 
F, that is 



T\ini...\ k ti k \ k+1 Hk+i S V [i] • • • V [i] F [* 



Each of these tensors describes the same physical properties of the electromagnetic field 
as the corresponding (super) fc -energy tensors do for the gravitational and massless scalar 
fields. They may arise at points of V n where the derivatives up to a given order of the elec- 
tromagnetic field vanish. From the discussion in sections || and || the dominant (supers- 
energy properties hold for the electromagnetic field, and one can define its (super) -energy 
densities and flux vectors relative to u. As in previous cases, we have 

Proposition 6.6 The (super) -energy (tensor) of the electromagnetic field vanishes at a 
point x £ V n if and only if the k-th covariant derivative of F is zero at x. In particular, 
all (super) -energy tensors of the electromagnetic field vanish in a domain D C V n if F 
vanishes in D. If the (super) -energy (tensor) of F vanishes in D, so do all the (super) - 
energy tensors with k > k. The manifolds containing a constant F are characterized by 
the vanishing of the super-energy (tensor) of F. ■ 
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In general (V n ,g), one can compute the divergences of the (super) -energy tensors 
to obtain very long formulae involving the Riemann tensor. However, in flat (V n ,g) it 
is obvious that Vm . . . Vm-Fra is completely symmetric in all the [l]-blocks, and thus all 
the corresponding (Ax^x)-pairs of the (super) fc -energy tensors can be interchanged in this 
case. This implies that, in a flat (V n , g), only two independent divergences of the (supers- 
energy tensors exist in principle. Computing them in this special case of flat manifolds 
one straightforwardly can prove using (p4l ) that in fact they are equal and vanishing, so 
that 



V7 rjpAlA*l...Afc/ifcAfc-|_ 1 /Zfc + l 



v. 



Vrn Fk 



in flat (V n ,g). 



Theorem 6.5 In any flat region of (V n ,g), all the basic ( super ) k -energy tensors of the 
electromagnetic field, as well as the general (k + l)(2fe + 1)!! -parameter family of tensors 
generated from it according to (jlqj, are divergence-free. ■ 



From this theorem one can define the following currents 



/ 



J\ 



V[i] . . . V[i] F[ 2 y, £2, ■ ■ ■ , 6fc+2 = T { 



(A/xi...A fc+1 /ife +1 ) 



Vrn • • • V m F 



[i] 



[1] *[2] 



tMl t-^fc+l fWs + 1 

S2 • • • ?2fc+l?2fc+2 



V 



which are divergence- free in flat manifolds (V n , g) for arbitrary Killing vectors £2* ■ ■ ■ > £2^+2 > 
leading to infinitely many conserved quantities in flat (V n ,g) for the electromagnetic field. 
This result is again not surprising, as explained for the massless scalar field, because if F 
satisfies the Maxwell equations (64) in flat (V n ,g), then the 2-forms Vi Ql . . . Vu,_ , , \F for 



\ a k-\ 



fixed values of a±, . . . , ctk+i are also solutions of the same equations, and their respective 
"energy-momentum" tensors of type (^) give rise to infinitely many conserved quantities. 

6.3 Massive fields 

The previous s-e tensors have been defined for the cases where the physical fields are 
massless. However, sometimes these fields carry mass, denoted generically by to, and 
then the field equations as well as the energy-momentum properties change. In this 
subsection the treatment of the massive cases is explained in detail and the corresponding 
(super) -energy tensors are derived. Only the prominent cases of the scalar and Proca 
fields will be considered explicitly, but the general rules will be valid for any other fields. 
In general, the existence of the mass amounts to having another quantity of physical 
relevance denoted by to which can be assumed to have (in natural units) dimensions of the 
inverse of length, L^ 1 . Thus, at first sight one realizes that the derivative of the field (say) 
has the same physical dimensions as the field multiplied by m, so that they are on the same 
footing as per how they can contribute to the energy or the higher order super-energies. 
Nevertheless, the derivative of the field has one more index block than the product of the 
field by m, so that one has to use the metric in order to combine them appropriately. 
Fortunately, the metric satisfies (up to sign) the DSEP, and thus one can construct new 
s-e tensors containing the mass which keep this fundamental property. As an illustrative 
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example, before giving the general definition, consider the energy-momentum tensor of a 
massive scalar field (ft, 

T Xfl = V A 0V^ - ^g Xi y p (ftVP(ft - ^g x „m 2 (ft 2 (70) 

which is symmetric and identically divergence-free if the Klein-Gordon equation 

V p V p (ft = m 2 (ft (71) 

holds. Let us observe that this tensor can be written as 

T Xfl = T Xfl { V[i]<£] + T {rri(ft} (-#A/J- 

The key idea here is that both Vm^ and rri(ft have the same physical dimensions, but of 
course the first is a 1-form and the second is a scalar. This mathematical difference is 
corrected by use of the fundamental metric tensor g. Other important points are: first, 
due to Property 1_1 the above tensor satisfies the DSEP because both 7\ M < Vm^} and 



T {m(ft} so do; and second, the divergence of the tensor Ta m is zero due to ( |7l| ) 

Bearing all this in mind, the generalization to arbitrary fields can be given in an 
inductive manner as follows. Assume that some physical massive field is described by 
an r-fold form £[ ni ],... ,[nr] (then, the 'smallest' s-e tensor of type (||) that can be formed 
will have 1r indices, so that this will be at the (super)' 7 '" 1 - 1 -energy level). Imagine then 
that we know the (super)' fe+r_1 - ) -energy tensor of this field which, by consideration of the 
physical dimensions, will involve basic s-e tensors of type (|6|) for i) the derivatives of the 
field of order k; ii) products of m by derivatives of the field of order k— 1; iii) products of 
m 2 by derivatives of the field of order k — 2; and so on. In other words, there will appear 
the s-e tensors of type (||) for a definite and precise set of tensors of type 

™*~%V],~.[«r] (72) 

where £ can take values £ = 0, . . . , k. Then, the basic (super) ( fc+r )-energy tensor of the 
massive t[ ni ],...,[n r ] is the one constructed adding all the s-e tensors (|]) for the derivatives 



V[l] (™ fc - / Vf L ]*[n 1 ],...,K]) 



of precisely the tensors appearing in the set (f72|), plus all the s-e tensors of type (||) for 
the products of m with exactly those in the list ( f72|) 



m [m k M,t M 



It is important to remark that not all possible tensors of the type m V fiir^i \ n j will 



necessarily appear in the list (72), and thus, in the next step, one only has to consider 
those arising from the ones really appearing in the list by taking the covariant derivative 
and multiplying by m. Observe further that the position of the m and of the derivative is 
important, and for instance the tensors wV[i]tr ni ] j ...r nr ,] and Vh](mtr„ 1 ] j .„ r^i) will produce 
s-e tensors with different index order. 

Some explicit examples will make all this clear. Starting with the scalar field (ft, at 
the energy level one can use the derivative of (ft as well as the product of (ft with m, or 
schematically 

V[i]0 m<ft, (73) 
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as tensors to construct the corresponding expressions (g). As we have seen in (|70|), this 
gives the correct energy- momentum tensor. Notice that in the first case the index in Vm 
gives rise to the pair of indices in the energy-momentum tensor, while in the second case 
it is m (so to speak) that gives rise to these indices by means of the corresponding —gx^- 
Now, in order to get the super-energy tensor of the massive scalar field, and following the 
general rule explained before, one must use the derivatives of the tensors in the list (|73[ ) 
plus the products of the mass m with the tensors in that list, that is 

V[i] (v [i]</> e m(j)j e m (y [i](/> © mA 

which produces the new explicit list 

V[i]V[i]<?f> V[i] (rrup) © rn (V[i]^>) ffimm</>. (74) 

Observe that the second and third terms are different because they produce different 
indices in the super-energy tensor. Explicitly, with the four tensors in the list (f74|), the 
following s-e tensors are built (all of them at the same level from the physical dimension 
viewpoint) 

Tap^ |V[i] V[i]</>| , T Q/3 |V[i] (m<f))j, T Xtl \m (v^Jj , T{mm<p} 

and the sum of all of them, including a —g for each m, produces the basic s-e tensor of 
the massive scalar field [Q: 

Sap\[i = T a pxn {V[i]V[i]</>| + T af3 |V[i] (m(j))j (-g x ^) + T Am \m (Vjij^Jj (-g a p) + 
+T{mm<t)} {-g a p){-9\n) = 2V a V {x cbV ^V ^ - g af3 (v A V^V M V p( /> + m 2 V x <f>V „cf>) 
-<?a m (V« V^V/jVp^ + m 2 V a 4>Vp(j>) + -g a(i g Xii (v^V^V^V^ + 2m 2 V^V'0 + m 4 2 ) (75) 



which has the same symmetry properties as (pOf). This tensor has been also found in |9 
for the case of Special Relativity (n = 4). 

Similarly, the (super) 2 -energy tensor for the massive (f) is obtained by using the list 
produced by taking the derivative and multiplying by m the tensors in (|74|), that is 



V 



[i] 



V[i]V[i]0© V[i] (m(t>) ffim (V[i]0J ®mm<j) 







© m V[i]V[i]0© V[i] (rruf)) ®m (V^J ®mra^ . 
Therefore, the basic (super) 2 -energy tensor of the massive scalar field is given by 

TafiX^rv \ V [1] V [1] V [1]^} + T af3X l i, { V [l] V [l] ( m( t>)} {-9rv) + 

+T a /3 TV [V[i] (mV[i]<^J| (-g\p) + T a/3 {V[i](mm<?!>)J (-gx^)(-g T u) + 

+r A/iTi/ [mV[i]V[i]^} (-5q/3) + ?V { mV [i](W>)} (-3a/3)(-5-n/) + 

+T Ti , |mmV[i]^| (-g a p)(-gxn) + T {mminc/)} (-g a)3 )(-g Xfl )(-g Tl ,) 

whose explicit expression can be read using (63), ( |60[ ) and (|57|). The (super) ^-energy 
tensors for the massive scalar field can be thus constructed for any k in this way. 
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Consider now a Proca field, also considerable as a massive electromagnetic field [10, 55]. 
This is a field satisfying the Proca field equations 

F pv = V M A„ - V u A p ( =► V [T F H = ), V P FP" = m 2 A v (76) 

from where one immediately deduces 

m ^0 =^ V U A U = 0. 

Therefore, in this case there is no gauge freedom and the 'potential' A is well-defined. In 
this sense, the tensor (^|) has a precise meaning, but unfortunately it does not have the 
proper physical dimensions comparable to those of (pq) . However, the existence of mass 
makes it possible to construct the corresponding tensor ( |65[ ) but using m A instead of A. 
This has the proper physical dimensions. Combining the two tensors we obtain 



T x ^ = T x ^{F [2] )+T Xp {mA {l] ) 



= F Xp F/ - jg^FpeFP" + m 2 A x A p - hm 2 g Xfl A p A p 

which is the standard symmetric energy-momentum tensor of the Proca field fll0|| . This 
tensor is divergence- free due to the Proca equations (|76|). 

Hence, the list of tensors used at the energy level has been found, namely 

F [2 ] © mA m 

and thereby, the list for the higher-order (super) fc -energy tensors can be deduced step by 
step. At the super-energy level, the list is given by 

V[i] (f [2] © mA[i|) m (F [2] mA [x ^j 

or equivalently 



[2] © Vji] [mA^A © mF[ 2 ] © mmA^. 



This example clarifies also that terms of type Fwm or mAmm (order is important because 
of the places of indices) will never appear according to our general rule. If one included 
those type of terms, the s-e tensor would not be divergence-free in flat (V n , g), see below. 
In summary, the basic s-e tensor of the Proca field becomes, by using (|67|), (|66|) and flBq) 



£af3\v = T a/3 Xfj, {V[i]F[2]| + T af5Xp |V[i] (mA^J } + T XfjL [mFp]] (-gap) + 
+T Xp [mm A {1] } (-g af3 ) = V a F Xp VpF/ + V a F^Vf3F x p - g a pV a F Xp \7 a F^ 



p 



-gxpVaFepVpF'"' + \g a pg Xl S 7 rF ap V ~F ap + m 2 [ V a A x V p A p + V«^ m V^a- 






9a^pA x V p A p - g Xp y a A p V p AP + -g a pg Xi y aA p V a A p 



-m 2 g aP ( F Xp F/ - ^g Xp ,F pa F p A - m 4 g a/3 U x A p - \g Xp A p A p ) . (77; 
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The next step is the (super) 2 -energy tensor of the Proca field, constructed using the list 
obtained by doing 



V[i] (V^-Fpj] V[i] [mA {1] j © mF [2] mmA {1] ^ 
© m (V[i]F[ 2 ] © V[i] (mAj © mF [2] mmij!]) 
which can be written as 

TapXfj.ru {V[l]V[l]F[2]| + T a/3 X/iTV |V[l]V[i] (mA[i]jj + 

+T a p TU |V[i] (mF [2 ]J J (~5A/i) + T a/3ru {V[u (mmA[i]J } (-SA/J - rn 2 g af3 £ XllT u 

with an explicit expression easily written by using (|69|), ( |77| ) and (@). All (super) ^-energy 
tensors are built similarly. 

Due to Property |4.1| , all the (super) ^-energy tensors constructed in this manner for the 
massive fields satisfy the DSEP, and obviously one can define the positive (super) -energy 
density, the corresponding flux vectors, etcetera. Moreover, the study of section || applies 
and one can construct, for instance, a 3-parameter family of s-e tensors for the massive 



scalar field starting from (75), and a 6-parameter family of s-e tensors for the Proca field 
starting from (JTJ\j. All this is obvious and will be skipped here. 

More important is the question of the divergence of the above tensors. For instance, 
for the s-e tensor (75) one gets exactly the same formula ( |dT| ) 

V a S> M = 2V /3 V (A 0i? M)p V^ - gx^R^VpV^V^ - 

- V ff </> (2 V V (A Bforf + gxpBfa W T 0) (78) 



and similarly for other cases. Thus, theorems ^^, |6.3| , |6.4| and 6.5 remain true for the 
massive fields. 

Theorem 6.6 In any flat region of (V n ,g), all the basic (super) k -energy tensors of the 
massive fields, as well as the general families of tensors generated from it according to 



(28), are divergence-free. 



This leads to the conserved currents in flat (V n ,g) using the same definitions as in the 
previous subsection. Details again are trivial and will be omitted here. 

Finally, there is a small change in the results proven in Propositions ^^ and |6.(3| . 
because the existence of mass implies that, at the (super) ^-energy level, all derivatives of 
the given field up to a given order (and not only those of this order) are involved, and 
thus the vanishing of the corresponding (super) ^-energy density implies the vanishing of 
all the derivatives up to the given order. More precisely 

Proposition 6.7 The ( super ) k -energy (tensor) of the massive scalar field vanishes at a 
point x &V n if and only if (p and all its derivatives up to the (k + l)-th order are zero at 
x. Analogously, the ( super ) k -energy (tensor) of the Proca field vanishes at x if and only 
if the k-th covariant derivative of F as well as A and all its derivatives up to the k-th 
order vanish atx. In particular, all ( super ) k -energy tensors of the massive scalar field, or 
of the Proca field, vanish in a domain D QV n if 4>, or respectively A, vanish in D. ■ 
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7 Exchange of super-energy quantities: conserved 
currents 

We finally arrive at the fundamental question of whether the s-e tensors have any phys- 
ical importance or describe some kind of relevant physical quantity. Obviously, the s-e 
tensors generalize the typical energy-momentum tensors, which are assumed to describe 
the energy-momentum and stress properties of a physical system. Hence, a naive guessing 
would lead to consider the s-e tensors as describing some kind of generalized energy- 
momenta or stresses of the same physical system. As far as mathematics are concerned, 
all the good properties of the energy-momentum tensors have been kept for the s-e ten- 
sors, namely, the positivity of the s-e density, the causality of the s-e flux vectors, and in 
general the DSEP which, as mentioned in section f|, allows to prove very general results 



on the causal propagation of fields and its uniqueness |1C, IS]. Let us also mention that 
the DSEP may serve to prove stability results such as those in [26 1 for the global stability 
of flat spacetime, but generalizing them to cases with matter or other asymptotic condi- 
tions. In general, the DSEP of the (super) ^-energy tensors provides estimates and a priori 
bounds on all the components of a given tensor and its derivatives of any order. Thus, 
from the mathematical point of view there is little doubt that the s-e construction herein 
presented is relevant. 

But, what about physics? Energy-momentum tensors are not only important from 
the mathematical point of view, but also from the physical one. The main role they play 
is in the definition of conserved quantities and its interchange between different physical 
systems. Thus, the main question arises: can the s-e tensors give rise to conserved cur- 
rents? And more importantly, can s-e quantities be exchanged between different physical 
systems? 

The first question has been already answered in this paper. As has been seen in 
section ||, the (super) fc -energy tensors of the scalar, electromagnetic or Proca fields provide 
infinitely many divergence-free currents in flat spacetime. Similarly, the generalized Bel 



tensor was proved to be divergence-free in Einstein spaces (R^ = Ag^), see Theorem 6.1. 
This last property leads to conserved currents if there exists at least a Killing vector, for 
due to the symmetries ( |3l| ) of the generalized Bel tensor only one independent divergence 
can be computed, and moreover, one has for arbitrary Killing vectors £i,£2>£3 



h 



[R[2],[2Y,£,1,{,2,£,3) = -B(a/3A) M C??2 ?3 = B (af3\v) Cf^2 £3 ( 79 ) 



so that from Theorem 5.1 it follows 



V M j M (-R[2],[2];£i,6,£3) =° in Einstein spaces: R^ u = Ag i 



[IV 



As already remarked after formula (|59|), not only the divergence- free property of B is 
needed, but also its symmetry properties which allow to form a divergence- free vector with 
any given Killing vector. It is noteworthy that sometimes ( |79| ) is divergence-free under 
much more general conditions (such as for instance when the Killings are hypersurface- 
orthogonal, and many others, independently of the form of the Ricci tensor), but these 
special cases fall out of the general treatment considered in this paper and will be presented 
elsewhere. 

Recapitulating, the generalized Bel tensor provides conserved currents if R^ = Ag^ u 
(i.e., in absence of matter), and the s-e tensors of the physical fields give conserved cur- 
rents in flat (V n , g) (that is, in the absence of gravitational field if this is described by the 
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curvature). This situation is, to say the least, encouraging, and one is led to the natural 
question of whether or not these respective currents can be combined to produce new con- 
served vectors involving several fields, such as gravitation and others. Let us remark that 
this situation is similar to that appearing in classical electrodynamics, where the energy- 
momentum tensor of the electromagnetic field is divergence-free in the absence of charges, 
but it is not if they are present. As is well-known, see for instance |30|, the divergence- 
free property (and a fortiori the conserved currents) is restored by adding the energy- 
momentum tensor of the existing charges, and the combined total energy-momentum 
tensor is the one without divergence, showing that the conserved quantities are the total 
energy-momentum and angular momentum of the charges plus the electromagnetic field. 
Alternatively, this can be seen as direct proof that part of the energy-momentum of the 
charges can be transferred to the field, and viceversa. This is what I call "exchange" in 
the second question asked above and the title of this section. 

Let us state, from the beginning, that the answer to this second question is positive, as 
is going to be shown explicitly in several cases. It is important to remark that, in order to 
find such mixed conserved currents, the corresponding field equations must be used. But, 
of course, the field equations will depend on the particular theory one is using (higher- 
order Lagrangians, scalar-tensor, string, dilaton, General Relativity,....), and thus the 
seeking of the conserved currents must be performed in each of these theories separately. 
In this paper, only the traditional case of Einstein's equations (in arbitrary dimension, 
units chosen such that 8ttG = c = 1) will be considered from now on. In particular, under 
this hypothesis, the interchange of super-energy between the scalar and the gravitational 
fields will be proved in full generality by constructing explicitly a divergence-free current 



from the sum of the generalized Bel tensor (pQ) and the s-e tensor (75) of the scalar field. 
Similarly, the exchange of super-energy between the electromagnetic and the gravitational 
fields will be shown by considering the propagation of discontinuities of Einstein-Maxwell 
fields. Nevertheless, before entering into this, let us make some remarks about other 
traditional treatments of the problem and the reasons of their failure. 

7.1 Non- validity of general divergence-free tensors 

The typical way to attack the problem has usually been the seeking of divergence-free ten- 
sors generalizing the Bel and Bel- Robinson tensors. Thus, soon after the publication of 
the Bel tensor, Sachs [86] found a four-index tensor T' a @ XfJ- in General Relativity (n = 4) 



which is identically divergence- free. This tensor was later generalized by Collinson (2£ 
who found the most general divergence-free tensor quadratic in the Riemann tensor and 
linear in its second derivatives if n = 4 . These results can be generalized to arbitrary 
n, see |35| (see also |J for related results). Unfortunately, these tensors do not satisfy 
any of the good mathematical properties of the s-e tensors @ herein presented, and their 
completely timelike component is not positive definite in general, nor they satisfy the 
DSEP. Furthermore, despite being divergence- free, all these tensors do not have any sym- 
metry property, and this invalidates them as candidates to produce conserved quantities. 
In order to have a conserved quantity one needs a divergence- free vector, so that Gauss' 
theorem is applicable. When one has a divergence- free tensor, in principle no conserved 
current can be built, unless there exists an auxiliary object which permits to find a true 
vector without divergence. Of course, this auxiliary object is usually a Killing vector (or 
some generalization thereof, such as conformal Killing vectors, Killing tensors, etcetera). 
The idea was explained in some detail when deriving (|59|). As remarked there, the sym- 
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metry property of the divergence-free tensor is crucial. Imagine, for instance, that we take 
a tensor such as that of Sachs (or Collinson), which is divergence- free 

V a T fal3Xfl = 

in the first index, but such that this first index has no symmetry properties with the other 
indices (in fact, Sachs' tensor is symmetric in the last three indices, T' a @ ^ = j 1/a (/ 3A M) ) 
but it has no symmetry involving a). Then, when seeking for a current j, and if there is 
at least one Killing vector £, one could try 

and then a direct computation gives, using that T" has no divergence and that £ satisfies 

dH) 

V Q j Q = 3T' a ^ (V^) £ A ^ = 3T'N^ (v^) £ A ^ + 

which is not zero in general, because the tensor T' a ^\ p is not symmetric in (a/5). The 
use of the completely symmetric part of T" (so that j a = r'( a/3A/i )£^£ A £ M ) is of no help 
here, because the divergences of T' a p\ p in any index other than the first is not zero. In 
summary, this type of tensors do not have the required mathematical properties which 
have been shown to hold in general for the s-e tensors (g) and, moreover, they may not 
be useful regarding the existence of conserved s-e quantities, which was the reason why 
they were constructed in the first place. 

Another example of this type is given by a four-index tensor presented in spinor form 
in [gQ] as the s-e tensor of the electromagentic field which, when combined with the 
Bel- Robinson tensor, was claimed to provide the appropriate s-e tensor for the case of 
Einstein-Maxwell spacetimes for n = 4, in General Relativity. The explicit expression of 
this tensor is 



P a ^ = 8 [V {p FP x V^F pa + V p F a{p V x F p f - 

-\g a {^xF pu V p) F ptJ + ^VaF pW V p F\g p)a - ^V p F aa V° F" \ p g Xp) 
which has the properties 

Papxn = Pa(^) , P af 5 P p = 0, V Q (V^ + P a ^) =0 in n = 4 

where in the last one the Einstein-Maxwell equations must be used and T is the Bel- 
Robinson tensor (|54|). Nevertheless, the tensor T + P has non- vanishing divergences in 
the last three indices, and there are no index symmetries involving the first index, so that 
it does not give divergence-free currents for general Killing vectors. Thus, for instance, 
P does not provide conserved currents for the rotational Killings in flat spacetime, in 
contrast to what happens with (^). The tensor P a p\ p is quadratic in the first derivatives 
V[i].F[2], but nonetheless it is not included in the six-parameter family E Q/ g AAt given in (|68|), 
not even allowing for negative constants ci,...,ce in (|68[), as can be checked. Actually, 
the tensor P a p\ p does not satisfy the DSEP nor the vanishing of its timelike component 
(which is non-negative) implies that Vm-Fpi is zero. All in all, it is doubtful that P be 
a good s-e tensor for the electromagnetic field and, in principle, the tensor ( |67[ ) (or its 
generalization E) is preferable. 

Now that we know the reasons why the general divergence-free tensors may not be 
good enough for the goal of s-e interchange, let us come back to the question of how to 
achieve this interchange and the conserved s-e quantities involving two fields. 

47 



7.2 Exchange of super-energy between scalar and gravita- 
tional fields: general divergence-free currents 

In this subsection an important theorem will be proved: if the Einstein-Klein-Gordon 
equations hold, then the sum of the s-e tensors (|3^) and ( |75| ) always provide conserved 
currents if there is a Killing vector M]. To that end, assume that there exists a (generally 
massive) scalar field 4> in the spacetieme (V n ,g) and that the Einstein field equations 

Rp, = V^V„0 + — ^mV V (80) 

hold, from where one can deduce the Klein-Gordon equation (|7l|). Then, the matter 
current appearing in (^) is 

2 

Jxfj.13 = V/3 V A V M - V/jV^ V\(p -\ —m 2 (j) (g/3f,Vx4> - gp\V ^4>) 

so that equation (37) becomes in this case 
V Q B a p Xfl = 2V^W (A <^R% )p/3 - g Xl >V*<f>V p V r <l>R°T P p + 2V (T V'V% )(t(a V m) </> - 

2V a <t>R a iXlM)f3 - 2V /3 0V A </>V M 0- 



n 



2 2, 

-m <ft 



-m 2 2 5/3(A V M) + g^Vtf \V p <jN p <t> + -mV 



n — I ' ' \ n 



(81) 

Similarly, expression ( |78D for the divergence of the basic super-energy tensor of the scalar 
field becomes now 

V a S a pXfi = -2V^VV (A 0i?% )p/3 + gx l y^V p V T <i>R\ P f 3 + 

+ fapWt + -l^mV) (2V^V (A V^ - g Xl NpV p (t) W) . (82) 



Therefore, the direct sum of tensors (30) and (|7q) is no£ divergence- free in general. But 
this is of no importance for our purposes, as will be proved in what follows. 

As remarked in the previous subsection, in order to find a conserved current one always 
needs an auxiliary object, usually given by a Killing vector, independently of whether the 
tensor providing the current is divergence-free or not. Thus, let us assume that £ is a 
Killing vector in our spacetime. A classical result [76, ^] ensures then that the massive 
scalar field is invariant under £, that is 

£"V M = O, (ifm^O). (83) 



If the scalar field is massless, then in fact one has £^V^</> =const. see [76, B2]. Whether 
or not the mass is zero, it always follows that 

e /3 V p (/>V /3 V p ^ = (84) 



where pq) and the derivative of (K3j) have been used. Contracting ( J8T| ) and (p2l) with any 
three Killing vectors £i,£2>£3 (or with three copies of the same if only one Killing vector 
exists) and using ((83-B5) one finds 



Pc^C^T? R a ,„ . - — Y7 a(^>T7 T7,.AU a P„ _!_„,. E> CT P„, T T7 Y7 A\ef 3 C X C^ 



£?£&£ V Q «S Vm) = - V a(p (2V p V (A J R CT / /3) + g^Rr'rfVpVrt) tf$% 
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and, in general, neither of these quantities is zero. However, as is obvious 

This still does not give a divergence-free vector, but now using again the symmetry prop- 
erties of (30) and ( f75| ) (which reveal themselves as essential in this calculation), one can 



rewrite the last formula as follows 







= £l £2^3 V " y B (.a/3\n) + S (ap\n)) = V " [ B (a(3Xfj,) + S (a/3Xfl)) £l ?2?3 

where (ESI) has been used. Therefore, the super-energy currents 

ja ( R [2],[2] <-> V[i] V[i]0; II, |*2, I3) = ( B (aP\») + 5 (a/3A M )) £1 £2 £3 ( 85 ) 

are divergence-free. Observe that only the completely symmetric part of the tensor is 
relevant here and thus, coming back to the discussion of section |5l perhaps only the 
symmetric part of the general tensor (28) is meaningful for physical applications. Notice 
that the vector ( |S5| ) reduces to the corresponding conserved current of the scalar field in 
flat spacetimes and to the vector (|79| ) if there is no scalar field. Thus, one can consider the 
vector field ( |85| ) as the appropriate generalization of those currents to the case involving 
both fields at the same time. 

Theorem 7.1 In any (V n ,g) satisfying the Einstein- Klein- Gordon equations $8(\ ) and 
containing at least a Killing vector, the current ma) is identically divergence-free. ■ 



As always, this produces conserved quantities by means of the Gauss theorem and the 
adequate integrals. The theorem also holds when m = for the case of the Killing vectors 



being either parallel or orthogonal to Vmi^. The relevance of Theorem 7.1 is that shows, 



explicitly, the interchange of super-energy properties between the gravitational and scalar 
fields, because neither Bt a p\n)£l ^2^3 nor ^(a/3X^)S,i^2^3 are divergence-free separately in 
the general mixed case, but each of them reduces to the corresponding conserved current 
in the absence of one of the two fields. 

7.3 Exchange of super-energy between electromagnetic and 
gravitational fields: propagation of discontinuities 

Special cases concerning the s-e interchange for the electromagnetic and gravitational 
fields have been considered in [61, ^, 98]. Due to the difficulty to handle the general 



case of an electromagnetic field with a Killing vector, mainly because the different possi- 
bilities arising (null or non-null electromagnetic field, with the Killing vector being also 



a symmetry of the electromagnetic field or not, etcetera [27, 37, 71, 73, ^lj) must be 
treated separately, in this paper only the illustrative and important case of the propaga- 
tion of discontinuities of the fields will be presented. This will be enough to prove the 
interchange of s-e quantities between the electromagnetic and the gravitational field. A 
thorough treatment of the propagation of discontinuities in general (including also other 
fields) will be given elsewhere [[|l]] (see also |8^, J9"0[), with the result that (super) ^-energy 
quantities and conservation laws arise naturally in this type of situations where the field 
has a 'wave-front'. 

For our present purposes, let us consider the case when there is an electromagnetic field 
K 2 i propagating in a background spacetime so that there is a necessarily null hypersurface 
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of discontinuity a [38, 62, |9l| (a 'characteristic'). Let us denote by [V] ff the discontinuity 
of any object V across a. Using the classical Hadamard theory [p2| , pq , 62] or equivalent 
procedures p8| , 91], one proves the existence of a 1-form c defined only on a and such 
that 

where n is a 1-form normal to a. Obviously, given that a is null, it follows that n is null 



n„n r 







so that n is in fact a vector tangent to a [pq| . Besides, n cannot be normalized, and thus 
it is defined up to a transformation of the form 



n 



yn, 



7> 



(86) 



where 7 is any positive function defined on a. The null curves tangent to n are necessarily 
null geodesies ('bicharacteristics') contained in a [38, ^, 91] 



n^V^n u 



Vn v . 



Obviously, an affine parameter could be chosen, using fl8q), so that \& can be set equal to 
zero for some n, but nevertheless the freedom (^6|) still remains partly. Observe that c is 
also affected by the freedom (|86|) as follows 

1 
c — > — c. 

7 

Let 5 denote the first fundamental form of <r, which is a degenerate metric in the null 
a because n is an eigenvector of g with zero eigenvalue, 53, 38, 87, 88]. Then, the second 
fundamental form of a relative to n can be defined by means of 



K. 



1 



i] 



£ ii 9ij 



where £^ denotes the Lie derivative [87] with respect to n within a. Kij is a 2-covariant 
symmetric tensor intrinsic to the null a and it shares the same degeneracy than g, that 
is, Kijn 1 = [37, S8|. Due to this property, even though there is no inverse metric 
for g, one can get the trace of the second fundamental form by contracting with the 
inverse of the metric defined on the quotient spaces defined by the tangent spaces to a 
modulus the degeneration vector n. This trace will be denoted by 1? and has the following 
interpretation: let S C a be any spacelike cut of a, that is, a spacelike (n — 2)-surface 
living in a and orthogonal to n. Then, 1? is a measure of the expansion of these cuts 
along the null generators of a as it can be easily related to the derivative along ft of the 
in — 2)-volume element of S [S8|. 

Making use of the Maxwell equations (^J) , the following propagation law (also called 
a 'transport equation' ||38||) is found (63, |63|, pO, 91] 



n^c 2 + c 2 (tf + 2^) = 0, c 2 



c^ > 0. 



This propagation equation ensures that if c^ = at any point x € a, then c^ = 
along the whole null geodesic originated at x and tangent to n. In fact, for arbitrary 
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conformal Killing vectors £1 , Q2 (or for two copies of the same if there is only one) , the 



(87) 



above propagation law allows to prove that |38| , 90. |91 

/ c 2 11^X1(2 w ls 
Js 



are conserved quantities along a, where u\g is the canonical volume element (n — 2)-form 
of S. What is meant by this is that 



2 /-^ /-V I 

c n^nXi (2 u\s 



2 rV- rv 1 



for any two spacelike cuts S and 5 of a. Of course, these conserved quantities can be 
easily related to the energy-momentum properties of the electromagnetic field because 
Tfj, u \ [F[2]]<t \ = c 2 n^n v , so that the energy-momentum tensor of the discontinuity [-F[2]] CT 

contracted with the conformal Killing vectors £i , (,2 gives the function integrated in (|87| ) . 
Further, this can be even related to the discontinuity of the electromagnetic energy- 
momentum tensor, because 



^{^[2]} = T^v \[F[2\]<t\ if tne background is empty. No- 



tice also that the conserved quantities (|87|) are invariant under the transformation (|8 



However, the integral (87) is trivial when [F t 



flU\ 



0, in which case c h 



0. When 



this happens, there arises a situation such that the electromagnetic field vanishes at a but 
its derivatives are non-zero there, so that the (super) -energy tensors must be relevant 



according to Proposition |6.6|. Using again Hadamard's classical results, now there exist a 
riant symmetri 
, |6l (6l p|, P, 



2-covariant symmetric tensor B pi/ and a 1-form f p defined only on a such that |j, |7[ flc| , 
4SL" 



[R. 



a/JA/iJo- 



B„ 



B v 



4ra [a- B /3][M n A]> "fw — ^vp, 

Again, these objects are affected by the freedom fl86|) according to 



2n»B^ + g^B^np 
0. 



0, 



n P f p 



(U,B, 



y.v ) 



2 \Ji^t>lw) 



Assuming that the Einstein-Maxwell equations hold, that is, (p4|) are valid and the Ein- 
stein field equations R pv — (R/2)g pv = T pu {F[ 2 ]} hold for the energy-momentum tensor 
( |6G| ) of the electromagnetic field, Lichnerowicz was able to prove the following propagation 
laws back in 1960 || (see also pfl) 

n^d^B 2 + B 2 ($ + 4* ) - 2 n a F ap B pT f T = 0, B 2 = B^B pu > 0, 
n^d^f 2 + f 2 ($ + 4tt) + 2n°F ap BP T f T = 0, f 2 = f^ > . 

Therefore, the discontinuities of the Riemann tensor and of Vm-Fra are mutual sources. 
This is yet another fact in support of our definition of (super) -energies and of the choice 
of levels: R]2], [2] an d Vm-Fra are at the same super-energy level, in this case the first k = 1. 
Actually, from d8l) and (p9h one clearly sees that Eol, p^, |63| 



n^ (B 2 + f 2 ) + (B 2 + f 2 ) (■& + 4* ) = 



so that, with the help of any conformal Killing vectors Ci> ■ ■ ■ > Ci ( or f° ur copies of the 
same one), the following conserved s-e quantities along a can be built 



B 2 + f 2 



P/-P, 



>Q >P /-H 1 ?-v i 

rianpn^rivCi C2 C3 C4 w |s 



(90) 
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in the sense that this integral is independent of the spacelike cut S. The fundamental les- 
son derived from this important formula is that both the electromagnetic and gravitational 
contributions are necessary, so that neither the integrals of B 2 n a n^n^n v C^Qi^z^i nor °f 
f 2 n a npn fJj n u Ci (2 (3C4 are equal for different cuts S in general. Another interesting point 



is that fl90|) can be related to super-energy tensors [|40|, |49|, |50|, |62j, |9l|| . The generalized 
Bel tensor ( |30| ) of the discontinuity [-R[2],[2]]o- becomes 

Tapiu, [[R[2},[2]\aj = 2B 2 n a n /3 n fl n 1/ 
while the s-e tensor ( |67| ) of the discontinuity [V[ 1 ]F[ 2 ]] CT reads 

T a/ 3^ u {[V[i]-F[2]] CT } = 2f 2 n a n p n^n u 
Therefore, the function integrated in (90) is simply 

2 [ T a0fiu {[-^Pl.pjlo-/ + T a/3)J.u { [V[i]-F[ 2 ]] CT }J C1C2 C3 Q 

which once again proves the interplay between super-energy quantities of different fields, 
in this case the electromagentic and gravitational ones. Observe that the above tensors 
are completely symmetric in this case, and that they together with the conserved quantity 
( |90| ) are invariant under the transformation (|86|), 

Actually, the exchange of higher-order (super) -energy quantities can also be demon- 
strated using this type of calculations. Again, the integral (p0) vanishes when f p = 0, 
B^ = 0, or equivalently when [V a F^ u ] a = [Rafix^L = 0. Thus, the (super) 2 -energy 
tensors are now relevant according to Proposition |6.6| . Under these hypotheses, one can 
prove that there exist a 2-covariant symmetric tensor A^ v and a 1-form Y p defined only 
on a such that plf 

[Vpi? Q/3 A M ] CT = 4,n p n [a A 0][fJi n x] , A^ = A VfJt 2rf A pfi + g^A^np = 0, (91) 

[V a V( 3 F llu ] a = 2n a npn [ ^Y u] , n p YP = 0. (92) 

These tensors are affected by (|S6|) according to 

{i/ii Apt,) ► — (i p , A^ v ) . 

Using again the Einstein-Maxwell equations one can prove transport equations analogous 
to (HP) 



n p d p A 2 + A 2 (i? + 6*) - 2 n a F ap A pT Y T = 0, A 2 = A pv A pv > 0, (93) 

n^d^Y 2 + Y 2 (0 + 6*) + 2n a F ap A pT Y T = 0, Y 2 = Y p Y p > . (94) 

Hence, the discontinuities of Vm-RpiM and of V^iV^F™ are sources of each other, which 
is again in favour of the definition of (super) -energies and of the choice of levels, in this 
case for k = 2. From (p3) and (E 



l d p [A 2 + Y 2 ) + [a 2 + Y 2 ) (0 + 6*) = 



so that, with the help of any conformal Killing vectors C,i, ■ ■ ■ , (e (which can be repeated), 
one finds mixed conserved (super) 2 -energy quantities along a similar to (|9(]) and which 
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need both the electromagnetic and gravitational contributions again. These new cur- 
rents can be related to the (super) 2 -energy tensors of the discontinuities [Vm-R^i^ltr an d 
[V[i]Vm-F[2]]er because, using ( |55| ) and ( |69| ) together with (91-^), one has 



Tapxuru { [V[i]-R[2],[2]]<t) = AA 2 n a n l3 n x nf l n T n u , 
Ta/jA^n/ {[V[i]V[i]F[2]] CT | = 4Y 2 n a n j3 nxn fl n T n u . 

In consequence, the integration over a of the functions 

^ \Taf3Xfiru {[V[l]i2[2],[2]]ff} + Tapxurv \\S [1]^ [l\ F [2\]a)) C1C2 C3C4C5C6 

provides the conserved (super) 2 -energy currents. Notice that this last expression, as well as 



the corresponding conserved quantity, is invariant under the transformation (86). Similar 
results can be derived for the general (super) fc -energy tensors of the discontinuities. 

8 Concluding remarks 



Very briefly, let us simply remark that the Definition 3.1 of basic s-e tensors has many 
good properties. From the mathematical point of view, they provide an appropriate 
way of defining scalars (the super-energy densities) and vectors (the s-e flux vectors) 
with the required mathematical properties, and especially they always have a definiteness 
property (the dominant property), which is strict for timelike vectors. Furthermore, the 
construction (y) is, on the one hand, universal (valid for any given tensor, in arbitrary 
dimensions, independent of any field equations), and thus all these good properties can 
be used with full generality; and on the other hand, it is fairly unique (apart from the 
index permutations studied in Section ||), as proved by the Lemma |4.2| and the Corollary 



4.3. The full uniqueness can be achieved by simply choosing the completely symmetric 
part, which in fact is the only one of relevance in the physical applications of Section g. 

Many mathematical applications are possible, such as i) the study of the causal 
propagation of fields and the uniqueness of the solutions, as noticed in Section |], see 



[H 16, 17, pig] ; ii) the use of s-e estimates to control the behaviour of the components 
of any given field and its derivatives ||, |l6| ; iii) the analysis of the non-linear stability 
of given spacetimes (V n ,g), or of any possible field therein contained, or both of them 
together. This may serve to improve and generalize the results of pq |; iv) the construc- 
tion of adimensional positive scalars (such as the qi,q2, <?3 of subsection |6.1|) measuring 



the quality of approximate solutions M, fl9|, or some physical feature as time-arrows or 



entropies [20, 79, 96 1, or the departure from some well-defined property of a given field; 



and many others. 



Of greater importance is that Definition 3.1 not only has mathematical relevance, but 



it is also physically very promising. The s-e tensors can be organized in a hierarchy of 
(super) fc -energy levels, but such that the tensors belonging to some level k for a physical 
field can be put in physical correspondance with the tensors of the same level for another 
field. Hence, the traditional super-energy Bel-Robinson tensor can be compared with 
the corresponding super-energy tensors of the electromagnetic or scalar fields. Thus, the 
correct physical dimensions for any given (super) -energy density seem to be those of 
energy density times L~ 2k . Furthermore, the analysis of the 'strength' of a field at points 
where its energy density vanishes but such that every neighbourhood of it contains the 
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field needs of the (super) -energy concepts, and this is why the s-e Bel-Robinson and Bel 
tensors arise naturally in General Relativity, where the energy density of the gravitational 
field can be always made to vanish at any point by appropriate choice of the reference 
system due to the equivalence principle. Analogously, the 'wave-fronts', shock waves and 
similar discontinuities must be analyzed from the (super) -energy viewpoint, specially 



those propagating causally along null hypersurfaces |4(], [49], g(], 62, ^33], 91, 99]. 

Finally, the most important point is that all the (super) -energy tensors give rise 
to divergence- free currents if the field generating the s-e tensors is isolated (that is, no 
other fields interact with it), and these currents can be combined to produce divergence- 
free currents mixing the s-e tensors of different fields, as explicitly shown in Section U 
in several outstanding cases. Therefore, the interchange of super-energy properties (in 
a wide sense: they can be super-momentum, or super-stresses etcetera) does happen, 
and the super-energy characteristics can be transferred from one field to another, such as 
energy properties do. A systematic study of the s-e exchanges between several fields in 
many different theories must be performed, and the explicit expressions of some of the 
conserved mixed s-e quantities should be analyzed in order to gain further insight into 
the possibilities opened by this new type of conservation physical laws. 
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Appendix 

Let (V n ,g) be any orientable n-dimensional differentiable manifold with a metric g of 
Lorentzian signature (—,+,. . • ,+)• The alternating Levi-Civita symbols are denoted and 
defined by 

X — X, . XMi-Atn _ slui-Hn] s _ r01...n-l _ , -i 

U /J,l...fl n — °[/il. ../*„]> ° — ° 1 "01.. .71-1 _ _ ' ' 

and then the canonical volume element n-form rj is given by 

Vltl...lin = %»!.../*»] = £ V l det Sl <W-Mn 

where e is a sign to be chosen according to the orientation. From here one can easily 
deduce 



„M1-Mn = „D*l~|*n] = c V det 9 

detg 
Therefore, one can write in general 



5" 



l—Hn 



„ yML.-Vn — I CLeX ff | <- vl .. Mn _ _sm...Vn (\ 1 \ 

where in the last expression the hyperbolic signature has been used and the so-called 
Kronecker tensor of order n has been defined 

SV\.„V n — r fy\—Vn 

U fil...fi n — U fJ.l-.-tJ.n U 
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In fact, one can define the Kronecker tensor of any order p < n by means of the following 
symbolic definition 



" Hi. -Up — 


XVl X"p 

"in ■ ■ ■ °Ml 
xvi s v p 

°H2 ■ ■ ■ °H2 




XVI SjVp 

°tlp ■■■ °flp 


so that 


xyi—v p _ r["l— "p] 
u lii-lip u ]pi...np] 


and one immediately gets 


c-Ul..M p U p+1 ...V p+q (P + <!)■ r.Vl...Up tVp+l — Vp+q \t 


n ... H , p ^ p+ i... fip+q ^ °lni...np°n p+ i...iip +q l 




cVl...VpfJ, p +l...lJ,p+ q 

0^ 1 ...^ p ^ p+ l...fM p+q 



p\q\ 



>+9J 

p+q ' 



A m 



lii... lip 



^m 



[ill... lip] 



p 






XVL.Mp a{(1] __ a{U} 

|Vl-Mj> m—v p — -n- in. ..fip 



(A.2) 
(A.3) 
(A.4) 



where A^- > lxx ___ l . denotes any tensor containing a set of p completely antisymmetric indices 
fix... fi p . 

With these ingredients the fundamental formulas (y-||) can be proved. As in the 
beginning of section g, let A^ ^ ll „, flpUl „ Mq = A^^!.,,^]^..,,,] denote any tensor with 
an arbitrary number of indices schematically denoted by {&}, plus two sets of p and q 
completely antisymmetric indices fi\ . . . fi p and v\...v q . Then, forming the duals with 
respect to these two sets of indices according to definition ([[]) 

/\{^} fq+l—Vn — „ „Vl~Vn A {U}tH—llp 

J ~ L * — I I ' fll-.-lln' ■"■ Vl...V„ 

llp + l-lln p\q\ ,pl Pn q 

and using here ( |A,1| ) and 



xyi—Vn J\\P.}lil--.lip 



pig! W-Mn 



VI... V q 



J^{Q} V q+ l...V n 

llp+l—fln 

U - sVl-Vq ?V g+ l...V n .{n}lll...llp 

p\(n - q)l(q\) 2 [Mi-Mg%+i-Mn] u i- u i 

from where, with the help of QA.4J) , one arrives at 



A {Q} ^ V q+ l. 

flp+l---lln 



n - j^mUl-Hv gV q+ l...V n 

p\( n — q)\n\ [lH-.-llq H q+ l...Hn} 



which is @. A similar calculation leads to @. 

Some particular cases of this identity are used in the main text. Let us start by 
considering the case in which p = q, and let us contract fi p +2 ■ ■ ■ fJ-n with v p +2 ■ ■ ■ v n , that 
is 

U - a {Q.} m... Up xWp+2-lln 

P)\(V\) 2 [Hl-Vp U \llp +2 -lln} 



/\{Q} llHp+2—tln 

^flp+2-.-lln 



(n 



which after use of ( |A.2| ) becomes 

^{O} HHp + 2..-lln _ _ 

<fyZp+2---£in 



7?! 



(n — p — l)l(p\) 2 
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A {Q}lll...llp ^[llrllp+2-llu] 

A [Ml-MP°A °iip +2 ...n n } 



and now, with the help of ( |A.3| ), this can be finally written as 

A{Q}pa 2 ...a p f A{^} W»p+2---Mn _ _XM /t{^}o-icr 2 ...(Tp 



(A.5) 
As a particular application of this general formula, we derive 

(p-1)! Aff2 -' Jp ' M + (n-p-l)! a Mp+ ;... m „,^ 



(P 



I /iM, O-2-..O-p | yl{n} cr 2 ...0-p 4{n}o-lCT 2 ...Cr p \ MM 

— 1V \ X(T2—ffp,fJ, ^ s± n Xcr 2 ...cr p iJXp-^- cr 1 o 2 ...a p I V A ' U / 



which has a righthand side symmetric in (Afi). Another application of (A.5) arises con- 
tracting A with fi, which leads to 

^{f2}o-io- 2 ...o-p _| j^{^} t p p+1 pp +2 ...p„ _ q (A. 7) 

p\ cr lC r 2 ...crp ( n _p)l Pp+iPp+2-P^n 



while contraction of A with fj, in ( A.6|) gives 

i 1 * 

/t{Q}aicr 2 ...ijp j /l{^} Pp+iP P +2---Pn 

(p-1)! aiff2 -' Jp + (n-p-1)! M, +1 /Vf2-Mn 



(2p — n) 



A {n}a 1 a 2 ...ap^ ^ 



I ~ ffitra-o-p 



so that the lefthand side of ( A. 6 ) is traceless in (A//) if n = 2p. 

Consider now any p-form S and let us apply (A.5) to the double (p,p)-form given by 
%n 1 ...n P %in...v P - One readily obtains 

1 1 1 

V, s^po 2 ...Op yi yifJ.fJ.p+2—^n A^V yi(Jl(T2...0-„ 

(A.9) 
Let us finally derive some formulas concerning any double (p, g)-form K^ 1 ..^ ptUl ... Uq = 

K \fn...lipUn...u q ]- From Q Jt follows 

.^ ^ J£0Pp+2---Pn,fJ-Cr2..-<T q _ 



(n — p — 1)1 ap p+2 ...p n ,Xa 2 ...a q 
1 1 

jrPp2—Pp,. IS IM?2...ffq Xfirs TS-plp 2 :.Pp,p,<T 2 :.<r q /A 1n \ 

/ _i\| XV \a 2 ...cr q I ^ap 2 ...pp, ^ u a 1 ^PlP2---Pp,Xo- 2 ...a q 1 ^- A^- iu / 

1 * 

_J( J^I3p 2 ...pp,pa q+2 ...a n _ 



(n — q — 1)! ap 2 ---Pp,Ao- 9+2 ...o-„ 

— K ^2—T q rc0P2—Pp,. xPzs jy-/3p2—p P ,(Ti-a- q fA ii\ 

— / _ -. w ■"■a/»2... ft,, " Ao- 2 ...o-«j ("A a P2— ft>i <T i— <T 8 A^- 11 / 

and double application of ( |A.5| ) gives also 

1 * * 

J£ J£PPp+2---Pn,p.cr q + 2 ---crn _ 

(n — q — l)!(ra — p — 1)! ap p + 2 ...p n Mq+2-<rn 

— I V-Pp2-Pp,P^2-ffqXr , XP'V;PP2--P V ,0\--<J q JS | _ 

~~ ( p _ 1)1(0 - 1)1 I J ^ap 2 ...p p ,X<T2-v q n* J ^ap 2 ...p p ,<j 1 ...a q I 

A/ 3 I fCPl-P P ,P^2-a- q ts XVjSpl---Pp,v\---v q K | fA 121 

Tjlffl _ 1)1 Q I -^pi...pp,ACT 2 ...cr 9 O^XV J^p 1 ... Pp ,CTl...ffq I ■ l^-^J 
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These three formulas ( [A.10 -A.12) are useful in finding the explicit expression of the super- 



energy tensor of Ky p ^y q y 

Let L a/ 3 : x^ = I'U/jijAd be any double (2, 2)-form (such as the Riemann or Weyl ten- 
sors). From (|3|) one has 

* ft ' 

V*m» _ 2!(n-2)!2! IAiw^...^] l A -"J 

and contracting here fi^ . . . /j, n with 1/5 . . . v n (if n > 5, of course), the following identity is 
found 

1 L , W..^» + X^ = i,^-^a^-^VS + ^Va ( A - 14 ) 



(n-4)! a/3p 5 ...p 
where 

-^aA — L a \ — — g a \ g L pa , L a \ = L a \ p . 

Eq. ( |A.14| ) is a generalization of the classical Lanczos identity |^8[ , and provides a direct 
relation between the double dual of Lw [2], ^[2], [2] itself and its traces. Notice that, if 
the double (2,2)-form Lr 2 i m is traceless, then the righthand side of ( |A.14 ) vanishes. The 



classical identity of Lanczos' [58] is for n = 4 and simply reads 



L.^+^^L^-L^-LVj+iVa. ifn = 4. (A.15) 
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